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ABSTRACT 


In  1979,  W.  B.  Colson  and  S.  K.  Ride  proposed  a  new  kind  of  electron  accelerator  using  a 
uniform  magnetic  field  in  combination  with  a  circularly-polarized  laser  field.  A  key  con¬ 
cept  is  to  couple  the  oscillating  electric  field  to  the  electron’s  motion  so  that  acceleration 
is  sustained.  This  dissertation  investigates  the  performance  of  the  proposed  laser  accelera¬ 
tor  using  modern  high  powered  lasers  and  magnetic  fields  that  are  significantly  improved 
over  the  values  in  the  original  analysis.  The  relativistic  equations  of  motion  for  the  elec¬ 
tron  are  solved  analytically  and  for  the  first  time  numerically  for  more  general  results  than 
previously  reported.  Realistic  beam  effects  are  investigated  to  determine  the  constraints 
on  initial  electron  energy,  injection  angle,  and  phase  to  achieve  optimal  performance.  This 
analysis  indicates  that  the  selection  of  an  electron’s  initial  conditions  are  not  as  stringent 
as  previously  thought.  Simulations  indicate  a  university  scale  system  can  achieve  electron 
energy  gain  of  over  one  GeV  in  a  meter. 
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CHAPTER  1 : 
Background 


In  1979,  W.  B.  Colson  and  S.  K.  Ride  [1]  proposed  a  new  kind  of  accelerator  using  a  uni¬ 
form,  magnetic  field  in  combination  with  a  circularly-polarized  laser  field.  Lasers  can  have 
strong  electric  fields  with  frequencies  high  enough  to  avoid  corona  formation  and  break¬ 
down.  The  key  is  to  couple  the  electric  field  to  the  electron  such  that  it  always  “sees”  a 
relatively  steady  field  so  that  acceleration  is  sustained.  Colson  and  Ride  [1]  suggested  that 
if  the  external  magnetic  field  is  aligned  along  the  light-propagation  axis  and  the  magnetic 
field  strength,  laser  frequency,  and  initial  electron  velocity  are  selected  judiciously,  the 
composite  field  forces  are  resonant  indefinitely,  so  that  the  electron  continues  to  be  acceler¬ 
ated  over  the  interaction  length.  For  laboratory  magnetic  fields  of  45  T  and  a  petawatt  laser, 
it  will  be  shown  through  numerical  simulations,  electrons  could  be  accelerated  from  about 
100  million-electron- volts  to  about  1  giga-electron- volts  (GeV)  in  about  a  meter.  This  is 
significantly  better  than  considered  in  the  original  research  paper  [1]. 

1.1  Motivation  and  Research  Questions 

In  today’s  high  powered  lasers,  electric  fields  on  the  order  of  10^^  V/m  are  possible  [2]. 
Electric  fields  are  transverse  to  propagation  and  oscillating  at  high  frequencies  on  the  order 
of  10^^  seconds  [2],  [3].  This  high  frequency  conspires  to  make  the  coupling  of  the  elec¬ 
trons  to  the  light  difficult.  Electromagnetic  forces  would  average  to  zero  due  to  the  rapidly 
oscillating  light  field.  Using  an  external  uniform,  static  magnetic  field  on  the  order  of  45 
T,  would  help  provide  sustained  coupling  between  the  electron  and  the  alternating  fields  of 
the  laser.  This  motivates  a  several  questions  that  this  dissertation  will  attempt  to  answer. 

1 .  Can  modern  lasers  and  magnetic  fields  provide  increased  coupling  compared  to  Col¬ 
son’s  and  Ride’s  original  results? 

2.  What  are  the  engineering  considerations  for  this  laser  accelerator  system? 

3.  How  is  success  measured? 

4.  What  are  the  constraints,  if  any,  for  optimal  performance? 
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5.  How  would  real  beam  effects,  such  as  emittance,  affect  performance? 


1.2  History  of  Particle  Accelerators 

Particle  accelerators  have  been  used  since  the  early  1900s  to  study  the  structure  of  matter. 
Since  then,  particle  accelerators  have  also  been  used  in  applied  research  for  industrial,  med¬ 
ical,  and  military  applications.  Particle  accelerators  are  broadly  classified  into  two  types: 
electrostatic  or  direct  current  (DC)  accelerators  and  oscillating  field  or  radio-frequency 
(RF)  accelerators. 

An  important  discovery  in  the  development  of  particle  accelerators  can  be  traced  back  to  the 
work  of  E.  Rutherford  [4].  In  1906,  Rutherford  developed  the  theory  of  atomic  scattering 
through  his  experiments  with  a-particles  and  atomic  disintegration  [4].  His  subsequent 
gold-foil  experiments  with  H.  Geiger  and  E.  Marsden  put  to  rest  the  “plum-pudding”  model 
of  the  atom  in  favor  of  the  nucleus  model.  The  Geiger-Marsden  experiments  proved  that 
unaccelerated  a-particles,  produced  by  the  spontaneous  decay  of  radium  atoms,  could  be 
scattered  by  over  90  degrees  when  fired  at  a  gold  foil.  Scattering  angles  of  this  magnitude 
could  only  occur  if  a  dense  nucleus  existed  at  each  gold  atom  site  [4]. 

By  1919,  Rutherford  believed  he  required  on  the  order  of  a  few  million  electron-volt  (MeV) 
ct-particles  to  split  the  lithium  nucleus.  This  value  exceeded  the  electrostatic  accelerator 
technology  at  that  time  so  research  in  this  effort  stalled  [5].  In  1928,  G.  Gamov  predicted, 
through  quantum  tunneling  effects,  that  approximately  450  keV  a-particles  would  suf¬ 
fice  in  splitting  the  lithium  nucleus.  Encouraged  by  Rutherford,  Cockroft  and  Walton  de¬ 
veloped  a  DC  accelerator,  shown  in  Eig.  1.1,  to  achieve  this  energy  requirement.  The 
Cockroft- Walton  accelerator  is  a  type  of  accelerator  whereby  static  electric  fields  are  used 
to  accelerate  charged  particles  [6].  Their  design  uses  a  voltage  multiplier  comprised  of 
capacitors  and  diodes  to  achieve  high  DC  voltage  from  a  low  DC  or  AC  voltage  source. 
By  1932,  their  accelerator  reached  400  keV  and  was  able  to  split  the  lithium  nucleus.  They 
were  awarded  the  Nobel  Prize  in  Physics  in  1951  for  this  work.  In  DC  or  electrostatic 
accelerators,  voltages  of  several  MV  can  be  reached  today  [5]  . 

Another  type  of  electrostatic  accelerator  is  the  Van  de  Graaff  generator  developed  by  R.  J. 
Van  De  Graaff  in  the  early  1930s,  see  Pig.  1.2.  In  this  design,  a  moving  belt  is  used  to 
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transfer  charge  to  a  hollow  metal  sphere.  Voltages  of  over  10  MV  can  be  reached  if  used 
with  an  insulating  gas.  Corona  discharge  limits  all  electrostatic  accelerators  to  voltages 
of  several  tens  of  MV  due  to  the  formation  of  ions  in  the  device.  The  stripped  electrons 
then  collide  and  cause  an  avalanche  of  more  electrons  leading  to  a  spark  discharge  in  the 
accelerator  and  thus  a  breakdown  of  the  electrostatic  field  [6],  [7]. 


Figure  1.1:  Cockroft-Walton  cascade  generator  used  to  split  lithium  atom,  from  [7]. 
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Conductor 


Figure  1.2:  Schematic  of  a  Van  de  Graaff  accelerator,  from  [7]. 


To  overcome  corona  formation,  G.  Ising  proposed  in  1924  the  use  of  radio-frequency  (RF) 
fields  [4].  Corona  formation  is  avoided  since  the  alternating  voltage  seen  by  molecules  is  of 
a  high  enough  frequency  that  the  electromagnetic  force  time  averages  to  zero.  In  1928,  R. 
Wilderoe  demonstrated  this  principle  with  50  keV  potassium  ions  [4].  Wilderoe’s  device  is 
a  linear  RF  accelerator  using  drift  tubes,  as  shown  in  Fig.  1.3.  The  RF  source  is  connected 
to  each  drift  tube  so  that  an  electron  only  feels  the  accelerating  force  of  the  electric  field 
at  each  drift  tube  gap.  Then,  while  the  electron  is  traveling  at  constant  velocity  in  the  drift 
tube,  the  fields  oscillate  such  that  the  electron  is  again  accelerated  in  the  next  gap.  Drift 
tubes  have  been  replaced  by  metal  cavities  since  the  mid- 1930s  and  by  superconducting 
RF  cavities  since  the  1960s  [4].  Superconducting  RF  cavities  minimize  the  heat  losses  in 
the  cavity  relative  to  the  energy  stored  in  the  resonator  achieving  a  very  high  Q-factor,  and 
therefore  high  fields  [4]. 
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Figure  1.3:  Schematic  of  a  linear  accelerator,  from  [7]. 


The  cyclotron,  developed  in  the  1930s  by  E.  O.  Lawrence,  makes  use  of  a  similar  accel¬ 
erating  structure,  but  is  more  compact  due  its  design  as  shown  in  Fig.  1.4.  Particles  orbit 
in  a  plane  between  two  poles  of  an  external  magnet  at  a  frequency  known  as  the  cyclotron 
frequency.  The  external  magnet  provides  a  static  magnetic  field  that  holds  the  particles  in 
the  plane  while  the  varying  RF  fields  on  each  of  the  D-shaped  electrodes,  called  “dees,” 
provide  the  accelerating  force.  The  charged  particles  accelerate  outward  from  the  center  of 
device  in  a  spiral  fashion.  Due  to  the  electron’s  small  mass,  the  cyclotron  was  limited  in 
energy  by  relativistic  effects  [4]  and  was  better  suited  for  heavier  particles  or  nuclei. 

A  special  type  of  cyclotron,  called  the  “synchrocyclotron,”  was  developed  by  E.  McMillian 
in  the  1940s.  In  this  design,  the  RF  fields  are  varied  to  compensate  for  the  relativistic 
effects.  Wilderoe  suggested  a  new  design  such  that  charged  particles  form  the  secondary 
loop  of  a  transformer  [4].  In  this  new  device,  called  a  “betatron,”  the  magnetic  field  is 
increased  as  the  particles  accelerate.  The  particle  trajectories  remain  circular  and  with  a 
constant  radius  [8].  Independently,  D.  Kerst  built  the  first  betatron  in  1940  that  produced 
2.2  MeV  electrons.  By  1950,  betatron  achieved  energies  of  300  MeV. 

To  achieve  higher  energies  and  to  maintain  a  constant  radius,  the  magnetic  field  had  to  be 
synchronously  adjusted  with  energy  of  the  particle  [7].  While  a  classical  cyclotron  uses 
both  a  constant  guiding  magnetic  field  and  a  constant-frequency  electromagnetic  field,  its 
successor,  the  “isochronous  cyclotron,”  works  by  local  variations  of  the  guiding  magnetic 
field,  compensating  for  the  increasing  relativistic  effects  of  the  particles  during  accelera¬ 
tion.  In  a  synchrotron,  this  adaptation  is  done  by  variation  of  the  magnetic  field  strength 
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in  time,  rather  than  in  space  [6].  In  proton  synchrotrons,  energies  of  over  1  TeV  can  be 
reached.  Some  energy  is  lost  through  synchrotron  radiation  that  is  emitted  when  charged 
particles  are  accelerated  radially. 

During  the  latter  half  of  the  last  century,  particle  accelerator  energies  nearly  doubled  every 
seven  years.  Since  most  exotic  particles  of  interest  do  not  exist  freely,  through  consequence 
of  Einstein’s  famous  equation  relating  energy  and  matter,  exotic  particles  must  be  created 
through  collisions  at  high-energies.  The  quest  for  higher  energies  to  unlock  new  exotic 
particles  has  pushed  accelerator  technology  over  the  last  century.  The  exponential  growth 
in  accelerator  energy  can  be  summarized  in  a  Livingston  chart  (see  Fig.  1.5).  Since  the 
1980s,  the  pace  of  accelerator  energies  has  fallen  off.  Today,  the  Large  Hadron  Collider  is 
the  most  powerful  accelerator  in  the  world.  Powered  by  a  proton  synchrotron  beam  line,  it 
can  reach  energies  of  seven  TeV. 

New  types  of  devices  using  plasma  waves  to  accelerate  electrons  are  being  developed. 
These  include  the  laser  wakefield  accelerator,  plasma  wakefield  accelerator,  plasma  beat- 
wave  accelerator,  and  self-modulated  laser  wakefield  accelerator.  These  technologies  may 
have  the  potential  to  accelerate  electrons  to  a  few  GeV  in  a  about  a  meter  [9] . 


Figure  1.4:  Schematic  of  a  cyclotron  accelerator,  from  [7]. 
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Figure  1.5:  Exponential  growth  of  accelerator  beam  energy,  from  [5] 
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CHAPTER  2: 
Laser  Accelerator 


This  chapter  will  discuss  the  proposed  components  for  a  laser  accelerator  system  using  real 
world  examples  to  frame  an  engineering  estimate  for  power,  size,  and  cost  estimates. 

2.1  Laser  Accelerator  System  Overview 

This  laser  accelerator  system  has  three  major  components: 

1.  Laser  source  -  provides  the  high  electric  field  used  in  accelerating  the  electrons 

2.  Electron  injector  -  produces  a  beam  of  electrons  and  directs  into  the  interaction  region 
at  an  initial  angle  and  energy 

3.  External  magnet  -  provides  the  magnetic  field  component  of  the  Eorentz  force  that 
enables  coupling  between  the  electrons  and  laser  light 

2.2  Laser  Source 

High  peak  power  lasers  with  large  electric  fields  can  be  laboratory  scale  and  provide  the 
force  for  significant  particle  acceleration.  As  an  example  of  high  peak  power  laser  tech¬ 
nology,  the  Texas  Petawatt  Easer  (TPE)  at  the  University  of  Texas  at  Austin’s  Center  for 
High  Energy  Density  Science  can  reach  energies  of  200  J  in  170  femto-seconds,  or  over 
1.1  petawatts  (PW)  peak  power  [2],  [10],  [11]  .  The  laser  system  provides  two  beam  lines 
directed  into  three  different  target  chambers  for  scientific  experiments.  The  main  beam  line 
is  based  on  optical  parametric  chirped  pulse  amplification  (OPCPA)  followed  by  power  am¬ 
plification  in  two  types  of  neodymium-doped  glass  [2].  Several  experiments  planned  for  the 
Texas  Petawatt  Easer  (TPE)  include:  production  of  fusion  neutrons,  study  of  laser  heated 
matter,  dynamics  of  shocked  materials,  radiative  hydrodynamics,  bright  x-ray  research,  and 
particle  acceleration  [10]. 

TPE  produces  170  fs  pulses  at  1057  nm  central  wavelength  [11]  .  Users  can  select  several 
energy  levels  from  60-200  J  per  pulse  at  a  rate  of  one  shot  per  hour,  or  7  J  per  pulse  at  four 
shots  per  hour  [11]  .  Table  2.1  summarizes  the  laser  system  specifications  [2]. 
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Table  2.1:  Petawatt  Laser  System  Speeifieations,  after  [2] 


Short  Pulse 

Eong  Pulse 

Pulse  Energy 

Pulse  Duration 

Wavelength 
Repetition  rate 
Eocused  intensity 

200  J 

170  fs 

1057  nm 

1  shot/hour 

>  10^1  W/cm2 

500  J 

2-20  ns 

527  nm 

1  shot/hour 

10^^  W/cm^ 

At  the  heart  of  the  petawatt  laser  system  are  three  OPCPA  stages  that  are  used  to  amplify 
the  seed  oseillator  energy  from  mJ  to  nearly  200  J  per  shot.  OPCPA  techniques  are  the  latest 
generation  of  chirp  pulse  amplification  (CPA)  technology.  CPA  has  enabled  scientists  to 
push  the  envelope  of  peak  power  systems  from  gigawatts  to  petawatts  in  the  past  30  years. 
CPA  was  first  demonstrated  by  D.  Strickland  and  G.  Mourou  in  1985  [3].  We  will  briefly 
describe  how  CPA  and  OPCPA  work. 

Consider  a  Gaussian  pulse  written  in  the  following  form  [3]: 

£(t)  (2.1) 


where  Eq  is  the  peak  electric  field  strength,  Up  is  a  measure  of  the  Gaussian  pulse  length,  bp 
is  the  linear  chirp  parameter,  and  (Oq  is  the  light’s  central  angular  frequency.  The  first  ex¬ 
ponential  factor  describes  the  Gaussian  envelope  and  the  second  term  describes  the  phase. 
The  total  instantaneous  phase  is  given  by  [3]: 


*^^total  —  -|-  bpt  . 


The  instantaneous  frequency  is  defined  as  [3] 


(Oiit) 


^0total 

dt 


Therefore,  the  instantaneous  frequency  for  the  pulse  given  by  Eq.  (2.1)  is  [3] 


(2.2) 


(2.3) 


(Oi{t)  =  C0o  +  2bpt. 


(2.4) 
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As  can  be  seen  from  Eq.  (2.4),  the  instantaneous  signal  frequency  increases  with  time  that 
is  called  an  up-chirp  for  bp  >  0.  Down-chirp  pulses  have  their  signal  frequency  decreasing 
with  time  with  bp  <  0.  Exponential  chirps  are  also  possible.  The  pulse  can  acquire  a  chirp 
during  transmission  through  a  dispersive  or  nonlinear  medium.  Eig.  2.1  shows  a  chirped 
Gaussian  pulse. 


Figure  2. 1 :  Diagram  illustrating  chirped  Gaussian  laser  pulse.  The  blue  curve  represents  the  highly 
chirped  signal.  The  red  curve  represents  the  Gaussian  envelope. 


Eig.  2.2  illustrates  basic  the  concept  of  chirped  pulse  amplification  [12].  An  oscillator 
produces  a  short  pulse  (1),  which  is  then  fed  into  gratings  that  stretch  the  pulse  temporally 
out  by  10^  -  10^  times  (2).  The  gratings  cause  the  low-frequency  components  to  travel  a 
shorter  path  than  the  high-frequency  components,  thereby  giving  it  an  up-chirp.  The  result 
is  a  temporally  stretched  pulse  with  a  low  enough  peak  power  to  safely  pass  through  the 
main  laser  gain  medium  usually  Nd:glass  or  Ti:sapphire  (3).  The  amplified  pulse  is  then 
recompressed  with  gratings  giving  it  a  down-chirp  resulting  in  a  very  short  and  high  energy 
pulse  (4). 
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1)  Oscillator  produces  short  pulses 


2)  Using  gratings,  the  pulse 
Is  stretched  by  10,000  x 


3)  Long  pulses 
are  now  sate  to 
amplify  In  laser 
amplifiers 


4)  After  amplification, 
the  pulses  are 
recompressed  with 
gratings 


5)  Resulting  pulse  Is 
short  (20 -1000  fs) 
and  high  energy  (>  1  J) 


Figure  2.2:  Diagram  illustrating  chirped  pulse  amplification,  from  [12]. 


Parametric  amplification  is  an  effect  where  the  dielectric  polarization  responds  nonlinearly 
to  the  electric  field  of  input  light.  A  higher  frequency  laser,  known  as  the  pump,  is  shined 
onto  the  parametric  crystal  resulting  in  a  splitting  of  the  pump  photon  into  a  pair  of  lower 
frequency  photons.  The  input  consists  of  a  pump  and  signal  photons  while  the  output 
consists  of  residual  pump,  signal,  and  idler  photons  respectively  such  that  COj  =  COrp  —  tOj. 
The  phase  relationship  between  signal  and  pump  determines  the  direction  of  energy  flow 
whether  there  is  amplification  or  deamplification  of  the  signal  [13] .  Several  benefits  of  para¬ 
metric  amplification  include:  a  much  wider  range  of  acceptance  wavelengths  compared  to 
laser  amplifiers,  higher  gain  per  unit  length  of  gain  material  in  pulsed  mode,  good  thermal 
properties  since  there  is  very  little  absorption  in  the  nonlinear  crystal.  Disadvantages  of 
parametric  amplification  include  no  pump  energy  acculumation  in  the  gain  medium,  pre¬ 
cise  pump  and  signal  synchronization,  and  losses  introduced  by  idler  photons  [14]. 

OPCPA  combines  the  techniques  of  CPA  and  parametric  amplification  to  achieve  very  high 
peak  power  in  ultrashort  pulses.  The  main  thrust  of  OPCPA  is  to  amplify  a  short  laser 
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pulse  while  mitigating  some  deleterious  nonlinear  effects,  such  as  self-focusing  and  phase 
modulation.  TPL  makes  extensive  use  of  OPCPA  technology  and  begins  with  a  100  fs 
tunable  oscillator  [2] .  The  pulse  is  stretched  to  2  ns  as  described  by  the  CPA  techniques 
above.  There  are  three  stages  of  parametric  amplification  using  nonlinear  crystals.  Stage 
one  consists  of  three  beta  barium-borate  (BBO)  crystals  and  a  Nd:YAG  pump  laser.  Stage 
two  uses  a  pair  of  BBO  crystals  also  pumped  by  a  Nd:YAG  laser.  Stage  three  uses  a  532 
nm,  4  J,  custom  laser  pump  and  two  potassium  deuterium  phosphate  (KDP)  crystals.  The 
final  preamplifier  stage  is  a  Ti:sapphire  crystal.  The  light  source  for  pumping  is  residual 
light  from  the  YAG  in  the  2nd  or  3rd  stage.  This  then  feeds  the  main  amplifier  comprised 
of  two  Nd:glass  amplifiers.  Pulse  energy  is  increased  to  about  300  J.  Finally,  light  is  passed 
to  the  compressor  where  it  is  housed  in  a  7x11  foot  vacuum  vessel.  A  multilayer  dielectric 
diffraction  grating  brings  the  pulse  energy  down  to  approximately  200  J,  150  fs  and  is 
directed  to  the  target  bay  for  experiments  [2]. 

OPCPA  enables  high  gain,  small  thermal  management  issues  (due  to  pulse  stretching), 
flexible  peak  gain,  and  spectral  gain  [10]. 

TPL  is  electrically  powered  with  a  large  capacitor  bank  next  to  the  laser/target  bay  rooms. 
Total  stored  energy  is  approximately  2  MJ.  Fig.  2.3  illustrates  the  layout  and  size  of  system 
with  the  laser  bay  and  target  area  comprising  of  approximately  1500  square  meters  [10]. 
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Figure  2.3:  A  schematic  diagram  that  illustrates  laser  bay,  target  bay,  and  capacitor  bank  rooms, 
from  [10] 


2.3  Electron  Injector 

Electron  injectors,  or  electron  guns,  come  in  many  different  varieties  and  configurations. 
The  basic  components  of  all  electron  injectors  are  the  cathode  for  emitting  electrons  into 
free  space,  a  beam  focusing  structure  usually  comprised  of  electromagnets,  and  the  subse¬ 
quent  accelerating  cavities.  Injectors  can  be  classified  by  their: 

1.  Emission  mechanism:  thermionic,  photocathode,  or  field  emission  cathode. 

2.  Electric  field  generation:  static  field  or  direct  current  (DC)  guns,  radio-frequency 
(RE)  guns,  normal  conducting  RE  (NCRE)  guns,  or  superconducting  RE  (SCRE) 
guns. 

Other  performance  properties  that  are  important  include:  average  beam  current,  peak  beam 
current,  electron  pulse  length,  emittance,  physical  size,  cooling,  reliability,  and  cost. 
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Thermionic  emission  of  electrons  uses  the  heat-induced  ejection  of  electrons  from  a  metal 
or  non-metal  surfaces.  A  common  example  of  this  is  a  vacuum  tube  where  a  hot  filament  is 
used  to  impart  sufficient  kinetic  energy  to  electrons  so  they  can  overcome  the  metal’s  work 
function  at  the  surface.  A  positive  anode  is  placed  some  distance  away  from  the  surface 
so  that  the  resulting  electrostatic  force  will  pull  the  electrons  into  vacuum  away  from  the 
cathode  where  they  can  be  further  accelerated. 

Photocathode  emitters  eject  electrons  from  the  cathode  by  the  photoelectric  effect.  A  drive 
laser  source  shines  light  energy  onto  a  metal  or  semiconductor  surface  [15].  The  quantum 
efficiency  (QE)  of  a  cathode  is  the  ratio  of  the  number  of  electrons  generated  to  the  number 
of  incident  photons.  Metals,  such  as  copper,  tend  to  have  a  low  QE  and  require  a  UV  drive 
laser  but  have  a  fast  response  time.  Semiconductors,  such  as  gallium  arsenide  or  cesium 
telluride,  have  higher  QE,  and  only  require  a  visible  drive  laser,  but  have  a  slower  response 
time.  The  UV  drive  laser  results  in  a  fast  response  time  while  the  visible  drive  laser  gives  a 
slower  response  time. 

Eield  emission  cathodes  eject  electrons  from  a  surface  due  to  a  nearby  high  electric  field 
gradient  that  lowers  the  work  function  so  that  electrons  can  tunnel  through  the  barrier  and 
escape.  Eield  emission  cathodes  are  also  considered  cold  cathodes  since  a  heating  element 
is  not  present.  Eield  emitters  structures  tend  to  have  spikes  or  sharp  tips  that  greatly  enhance 
the  electric  field  at  the  tip  near  the  surface.  Recent  papers  have  shown  that  carbon  nanotubes 
may  be  excellent  field  emitters  [16]. 

Eollowing  the  cathode  emission  section  of  the  electron  gun,  DC  or  RE  fields  provide  the 
initial  stage  of  acceleration.  RE  cavity  shapes  can  vary  from  pillbox,  elliptical,  quarter- 
wave,  or  spoked.  Spoke  cavities,  see  Eig.  2.4,  are  a  coaxial  half-wave  length  cavity  with  an 
outer  conductor  turned  ninety  degrees  so  that  its  axes  is  directed  along  the  beam  path  [17]. 

Electron  guns  with  a  photocathode  can  be  combined  with  a  DC,  NCRE,  or  SCRE  acceler¬ 
ating  cavity  section.  NCRE  gun  operations  span  a  broad  range  of  RE  frequencies  ranging 
from  a  few  MHz  to  several  GHz.  The  RE  field  and  the  drive  laser  pulses  must  be  care¬ 
fully  synchronized  so  that  the  electrons  ejected  via  the  photoelectric  effect  are  properly 
accelerated. 
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Figure  2.4:  Cross-section  of  a  triple  spoke  cavity,  from  [17] 


In  NCRF,  resistive  losses  eause  the  eavity  walls  to  heat  up  and  require  eooling  systems 
based  on  a  water  heat-exehanger.  Heating  is  greatly  redueed  through  use  of  a  supereon- 
ducting  eavity  but  system  eomplexity  and  cost  is  increased.  The  liquid-He  requirement  for 
some  super-conducting  cavities  may  be  somewhat  mitigated  by  choosing  a  lower  RF  band. 
For  RF  of  less  than  about  500  MHz,  a  liquid-He  system  of  only  4K  is  sufficient.  A  2K 
liquid-He  system  is  required  for  RF  greater  than  about  500  MHz. 

Electron  guns  can  be  continuous  wave  or  pulsed.  A  particle  synchronized  with  the  RF  field 
is  called  a  synchronous  particle.  Particles  bunched  near  the  synchronous  particle  in  the 
oscillation  of  the  RF  field  is  called  a  bucket.  Electrons  can  be  present  in  each  RE  bucket, 
or  in  pulsed  guns  where  certain  RE  buckets  are  filled.  Typical  frequencies  are  a  few  MHz 
to  several  gigahertz  with  bucket  charges  from  pC  to  nC  per  bunch. 

As  mentioned  earlier,  electron  beam  emittance  is  an  important  quality  to  consider.  It  is 
a  measure  of  particle  spread  in  position  and  angle  and  is  typically  measured  in  units  of 
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mm-mrad.  Normalized  emittance,  £„  is  defined  as 


£n  =  ro<r><  e  >  (2.5) 

where  jo  is  the  beam  Lorentz  faetor,  <  r  >  is  the  average  spread  in  partiele  position  and 
<0>  is  the  average  spread  in  partiele  angle.  As  we  will  see  later,  the  eleetron  beam  positions 
must  be  well  inside  the  optieal  mode  of  the  laser  to  be  aeeelerated.  For  a  relatively  good 
beam  emittanee  of  10  mm-mrad  and  Yo  =  130,  the  spread  in  position  and  angles  must  be 
less  than  0.08  mm-mrad.  We  shall  see  later  that  this  is  a  suffieiently  small  spread  in  angles 
in  order  to  aehieve  optimal  energy  gain. 

The  Naval  Postgraduate  Sehool  and  Niowave,  Ine  have  designed  and  developed  the  Mark  I 
SCRF  injeetor  (see  Fig.  2.5).  This  design  uses  a  quarter  wave  structure,  which  are  known 
to  achieve  a  high  accelerating  efficiency.  Furthermore,  operating  at  a  low  RF  band  of 
350  MHz  allows  use  of  a  4  K  helium  refrigeration  system  thereby  greatly  minimizing  the 
complexity  of  total  system.  The  Mark  I  uses  an  RF-coupler  with  a  nominal  electron  beam 
energy  of  a  few  MeV. 


Figure  2.5:  Naval  Postgraduate  School  SCRF  Gun 


17 


2.4  External  Magnetic  Field 

The  external  magnetic  field  in  the  laser  accelerator  system  serves  to  provide  coupling  be¬ 
tween  the  electron  beam  trajectory  and  the  laser’s  transverse  electric  field.  It  will  be  shown 
that  magnetic  field  strengths  on  the  order  of  30  T  to  60  T  provide  for  sufficient  electron 
coupling.  We  will  find  that  higher  magnetic  fields  allow  the  initial  electron  energy  and 
electron  injection  angle  to  be  much  lower. 

An  example  of  high  magnetic  field  technology,  consider  the  National  High  Magnetic  Field 
Laboratory  (NHMFL).  NHMFL  currently  holds  the  record  for  the  strongest  continuous 
static  field  in  the  world  at  45  T.  NHMFL  is  a  federal-state  research  facility  in  partnership 
with  Florida  State  University,  the  University  of  Florida,  and  Los  Alamos  National  Labo¬ 
ratory  under  a  National  Science  Foundation  grant.  Its  mission:  “To  provide  the  highest 
magnetic  fields  and  necessary  services  for  scientific  research  conducted  by  users  from  a 
wide  range  of  disciplines,  including  physics,  chemistry,  materials  science,  engineering,  bi¬ 
ology  and  geology”  [20].  There  are  several  magnet  types  available  for  scientists  to  conduct 
basic  research.  Table  2.2  summarizes  some  of  NHMFL’s  magnets  that  are  available  [18]. 

Table  2.2:  Available  Magnets  at  the  National  High  Magnetie  Field  Laboratory,  after  [18] 


Field  [T] 

Bore  [mm] 

Power  [MW] 

Resistive  and  Hybrid  Magnets 

20 

195 

20 

24.5 

32 

14 

25 

52 

19 

27 

32 

14 

30 

32 

16 

33 

32 

16 

45 

32 

15 

Superconducting  Magnets 

Temperature 

15 

45 

10  mK  -  1  K 

17.5/19.5 

52 

0.4  -  300  K 

18/20 

52 

20  mK  -  2  K 

The  Hybrid  Magnet  at  the  NHMFL  is  the  world  record  holder  for  continuous  magnetic 
field  (see  Fig.  2.6).  It  consists  of  a  room  temperature  magnet  (31  T)  surrounded  by  a 
superconducting  magnet  (14  T)  to  produce  a  continuous,  static,  45  T  field  in  its  32  mm 
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bore  over  one  m  length.  The  power  supply  is  rated  for  40  MW  [20].  The  eooling  system 
is  buffered  by  a  3800  eold  water  storage  tank  that  is  maintained  at  7°C  by  a  28  MW 
ehiller  [20],  [21]. 


Hybrid  Magnet 
Cross  Section 


Electrical  and 
Helium  Supply  for 
Superconducting 
Magnet 


Outsert  Magnet 
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Figure  2.6:  Cutaway  diagram  of  NHMFL  45T  Hybrid  magnet,  from  [21]. 


The  supercondueting  “outsert”  magnet  wraps  around  the  room  temperature  magnet.  It 
is  eomprised  of  three  eonneeted  subeoils  in  series  with  a  design  eurrent  of  10  kA.  Eaeh 
subcoil  is  made  with  of  a  niobium- tin  alloy  with  a  combined  stored  energy  of  nearly  100 
MJ.  Many  of  the  supported  research  areas  for  the  hybrid  magnet  include:  condensed  matter 
physics,  quantum  hall  effects,  high  temperature  superconductivity,  geochemistry,  optical 
microscopy,  resonant  ultrasound  spectroscopy  and  NMR  studies,  for  example.  The  system 
was  design  for  maximum  user-friendliness,  easy  access  to  bore  experiment  area,  and  a 
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10-year  life  with  2000  charge/discharge  cycles  for  the  “insert”  or  inner  restive  magnet  [22]. 


For  use  within  the  laser  accelerator  system,  static  and  continuous  magnetic  fields  may  not 
be  needed  because  the  time-scales  of  interaction  are  very  short,  typically  <  10  ns.  Pulsed- 
field  systems  are  able  to  provide  long  duration  pulses  compared  to  the  time-scales  of  inter¬ 
action  for  the  laser  accelerator  are  in  effect  “static."  NHMFL  Pulsed-Field  Facility  (PFF) 
at  Los  Alamos  National  Laboratory  (LANL)  offers  researchers  access  to  various  capacitor 
driven  magnets  up  to  fields  of  100  T  (see  Table  2.3)  [19].  The  fields  can  be  maintained  over 
a  length  of  about  one  meter. 

The  PFF  system  (see  cutaway  Fig.  2.7)  combines  a  40T  outer  magnet  with  a  60T  field 
created  by  a  capacitor-bank  insert.  The  power  system  consists  of  a  1.43  GW  generator 
combined  with  seven  64  MW  power  supply  modules  [23]  driving  each  of  the  magnet’s 
seven  coils. 

NHMFL  Pulsed  Field  Facility  is  currently  designing  the  next  generation  of  non-destructive 
pulsed  field  magnets.  The  lOOT  magnet  is  a  joint  project  between  the  National  Science 
Foundation  and  the  Department  of  Energy.  Pulses  up  to  lOOT  for  periods  of  milliseconds 
will  be  available  for  users.  Table  2.4  recaps  the  systems  we  have  used  in  our  examples. 

Table  2.3:  Available  Pulsed  Magnets  at  PFF  LANL.  SP  =  Short  Pulse.  MP  =  Mid- Pulse,  after  [19] 


Cell  No. 

Magnet 

Pulse  Duration  (ms) 

Bore  (mm) 

1 

50TSP 

25 

24 

2 

40TMP 

400 

24 

2 

50TMP 

400 

15 

3 

60TSP 

25 

15 

4 

65TSP 

25 

15 

5 

60TSP 

25 

15 

Table  2.4:  Magnet  Speeifieation  Summary,  after  [19]. 


System 

Power 

Cooling  Power 

Volume 

45T  Hybrid 
Pulsed 

25  MW  (continuous) 
1.43  GW  @  85  ms 

9  MW 

5  MW 

175  m^ 
35  m^ 
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Figure  2.7:  Cutaway  diagram  of  NHMFL  lOOT  pulsed  magnet  system,  after  [19] 
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CHAPTER  3: 
Equations  of  Motion 


3.1  Relativistic  Lorentz  Force 


We  begin  the  laser  accelerator  analysis  by  considering  whether  a  classical  or  quantum  me¬ 
chanical  treatment  is  required.  The  injector  gives  the  electrons  an  energy  of  about  120  MeV 
with  a  one  per-cent  energy  spread  as  input  to  the  laser  accelerator.  This  corresponds  to  a 
spread  in  electron  wavenumbers  of 

=  m-\  (3.1) 

he 

So  the  corresponding  uncertainty  in  position  is 

m.  (3.2) 

2Ak 

The  uncertainty  in  position  is  much  smaller  than  the  laser  wavelength  of  1  jim.  The  laser 
accelerator  is  clearly  in  a  classical  regime  per  Ehrenfest’s  theorem  [24]. 

Continuing  the  analysis  with  the  relativistic  Lorentz  force  law  in  the  general  case,  an  elec¬ 
tron  moving  through  a  magnetic  field  will  spiral  around  the  field  lines.  Without  an  applied 
electric  field,  the  electron  will  not  gain  energy.  The  strong  electric  fields  found  in  today’s 
high  powered  lasers  can  provide  huge  accelerating  gradients.  If  a  laser’s  electric  field  can 
be  coupled  to  an  electron’s  spiral  motion,  it  can  lead  to  significant  energy  gain.  In  compact 
covariant  notation,  the  relativistic  Lorentz  force  law  [25]  in  SI  units  is 


dp^  dU^ 

=m—— 
dT  dx 


-eF^^Uy 


(3.3) 


where  =  (ymc,  p)  is  the  four-momentum,  7=1/  a/I  —  (v/c)^  is  the  Lorentz  factor,  m  is 
the  mass  of  the  electron,  =  7(c,v)  is  the  covariant  form  of  the  electron’s  four- velocity, 
X  is  proper  time,  e  is  the  magnitude  of  the  electron’s  charge,  and  is  the  field- strength 
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tensor, 


f  0  -E^/c 
E./c  0 
Ey/c 

\  E^/ C  ~Ey 


Ey/c  E^/c  ^ 

—  B^  By 

0  -B, 

Bx  0  y 


(3.4) 


We  can  transform  between  covariant  and  contravariant  forms  of  4- vectors  in  the  usual  way 
using  the  metric  tensor 


S^v 


( 1 

0 

0 

0  ^ 

0 

-1 

0 

0 

0 

0 

-1 

0 

0 

0 

-V 

(3.5) 


The  left-hand  side  of  Eq.  (3.3)  is  a  proper  time  derivative.  In  order  to  transform  this  into 
a  lab-frame  time  derivative,  we  make  use  of  the  invariance  of  the  spacetime  interval.  For 
some  inertial  frame  K,  the  spacetime  differential  interval  is  given  by 


ds^  =  c^dt^  —  dx^. 


(3.6) 


For  another  inertial  frame  K',  the  interval  is 


ds^  =  c^dt'^  —  dx!^ . 


(3.7) 


For  sufficiently  small  spatial  intervals. 


dx  =  y{t)dt. 


(3.8) 


Substituting  Eq.  (3.8)  into  (3.6)  yields 

ds^  =  c^dt^{l-l5^)  (3.9) 

where  =  (/3x,  /3z)  and  =  v^/c^.  Since  V  <  c  in  all  situations,  ds^  is  always  greater 

than  zero  in  Eq.  (3.9).  For  the  given  metric  in  Eq.  (3.5),  positive  spacetime  intervals  are 
called  timelike.  Timelike  intervals  are  ones  where  there  exists  a  proper  time  frame  K'  such 
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that  the  spatial  interval,  dfd ,  is  zero.  So  we  have  from  Eq.  (3.7), 


ds’^  =  c^dt'^  =  c^dx^. 


(3.10) 


The  spacetime  interval  is  an  invariant.  Therefore,  Eqs.  (3.6)  and  (3.10)  are  equivalent. 
Combining  the  right  hand  sides  of  Eqs.  (3.9)  and  (3.10)  we  now  have, 

c^dx^  =  c^dt^{l-p^).  (3.11) 


We  can  now  perform  a  change  of  variables  in  Eq.  (3.3)  using  dt/dx  =  7  from  Eq.  (3.11), 


d  dt  d  d 

dx  dx  dt  dt 


(3.12) 


Therefore,  Eq.  (3.3)  becomes 

^  =  =  (3,13) 

dx  dt 

We  will  solve  Eq.  (3.13)  subject  to  the  laser  field  and  an  external  uniform  magnetic  field 
where  the  appropriate  forms  for  the  electric  and  magnetic  fields  will  be  added. 


3.2  Equations  of  Motion 

Consider  a  circularly-polarized  plane  wave  given  by  the  following  expressions  for  the  laser 
electric  and  magnetic  fields,  and  Br,  with  an  external  uniform  magnetic  field.  By,  aligned 
along  laser  propagation  axis,z. 


with 


Eo(cosC,-sin^,0), 

(3.14) 

—  (sin^,cos^,0), 
c 

(3.15) 

(0,0, 5„), 

(3.16) 

^  —  kfZ  cOft . 

(3.17) 
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The  plane- wave  phase  is  given  by  the  angular  frequency  of  the  laser  light  is  given  by  cOr, 
and  the  laser  light  wavenumber  is  given  by  kr.  Substituting  Eqs.  (3.15)  into  (3.4)  yields 


/  0 

-^cosC 

^sinC 

0  \ 

f^cosC 

0 

-Bu 

^cosC 

-fsinC 

Bu 

0 

-fsinC 

V  0 

-^cosC 

^sinC 

0  / 

(3.18) 


for  the  field-stress  tensor.  Substituting  Eq.  (3.18)  into  (3.3)  and  running  through  compo¬ 
nents  /i  =  0, 1, 2, 3  results  in  the  following  equations: 

^  =  -—[pxCosC-pySinQ  (3.19) 

at  me 

4(rfc)  =  -  — (3.20) 
at  me  m 

=  — (l-i8,)sinC  +  ^i8.  (3.21) 

at  me  m 

d  pF 

;7:(r^z)  = - -  [I5x  cos  C-jiy  sin  Q  (3.22) 

dt  me 

where  f^x  —  Vxjc,  Py  =  Vy/e,  and  /3^  =  v^/c.  Eet  us  describe  the  geometry  carefully  (see  Eig. 
3.1).  Given  initial  electron  energy  we  have  the  following  definitions: 


i8o  = 
/3zo  = 
Pxo  = 
Pyo  = 


l3o  cos  a  = 
[5  cos  (j)o 
15^0  sin  (jto 


1 

To 


irx/yl -^^osa, 


j5o  sin  a  cos  (jto 
Po  sin  a  sin  (j)o  -- 


1 

Yo 

1 

—  1 
Yo 


yI- 


—  1  sin  a  cos  (jto,  and 


1  sin  a  sin  (j)o 


(3.23) 

(3.24) 

(3.25) 

(3.26) 


where  a  is  injection  angle  of  the  electron  measured  from  the  z  axis  (optical  and  magnetic 
field  axis),  j8o  is  the  initial  electron  speed  relative  to  the  speed  of  light,  is  the  x  compo- 
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nent  of  the  initial  electron  speed  relative  to  the  speed  of  light  (similarly  for  /3yo  and 

is  the  perpendicular  component  of  the  initial  electron  velocity  relative  to  the  speed  of 
light,  and  (j)^  is  the  initial  angle  between  /3^o  and  /3xo- 


Z 


laser  light 


Figure  3.1:  Components  of  the  laser  accelerator  interaction.  The  electric  field  vector  is  given  by 
Ef.  Velocity  vector  (in  units  of  c)  is  described  by  )3  with  components  Px’  Py’  and 
Transverse  velocity  is  given  by  j3j^.  The  injection  angle  is  given  by  a.  The  angle 
between  the  rotating  electric  field  and  P^  is  defined  as  6.  The  initial  azimuthal  angle, 
00,  is  defined  as  the  angle  between  the  x  axis  and  at  t  =  0.  The  external  external 
magnetic  field  is  given  by  Bu. 


Notice  that  the  right-hand  side  of  Eqs.  (3.19)  and  (3.22)  are  equivalent.  Equating  the 
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left-hand  sides  to  each  other  and  integrating  gives  a  constant  of  motion, 


dy 

dt 


r(i-^z)  =  e 


(3.27) 


where  e  is  a  constant  for  all  time  [1].  Both  yand  are  dynamical  variables.  As  yincreases, 
the  factor  1  —  must  decrease  to  keep  e  constant.  We  can  express  -h  in  terms 

of  e  as  follows.  From  the  definition  of  the  Lorentz  factor  we  have, 


7  = 


(3.28) 


Using  Eq.  (3.27)  for  results  in 


I3l  =  ^l2re-e^-l].  (3.29) 

From  Eqs.  (3.27)  and  (3.29),  as  the  electron  gains  energy  from  the  laser  accelerator,  y 
increases.  Therefore,  /3^  will  increase  while  /3^  decreases.  The  coupling  to  electric  field  is 
governed  by  . 


Now  we  will  attempt  to  uncover  other  constants  of  the  motion  by  extracting  a  full  time 
derivative.  Consider  the  following  relations. 


d 

dt 


sin^ 


—  sin(k^z—  (Ort) 

COs(  krZ  —  COrt)  —  {Kz  —  (Oyt) 
dt 

COs{krZ  —  COrt)  {Kz  —  (Or) 

—  (Or{l  —  Pz)  COs{krZ  —  COrt) 
-COr{l-  I5z)  COS  C 


(3.30) 


and  similarly. 


d 

dt 


cos^  =  (Or{l  —  j8z)sin^. 


(3.31) 
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Furthermore,  applying  the  chain  rule  on  J/Jt(7/3xSin^)  and  substituting  Eq.  (3.30)  yields 


^  irl^x  sin  C)  =  sin  C  j^iYPx)  +  C) 


=  sin  C  ^  (r/3x)  -  r^x®r(  1  -  ^z)  cos  c , 


(3.32) 


and  also  using  Eq.  (3.31), 


^  (rPy  cos  C)  =  cos  C  ^  (rPy)  +  rPyG)r(l  -  ^z)  sin  C- 


(3.33) 


A  second  constant  of  motion  can  be  determined  by  solving  Eqs.  (3.20)  and  (3.21)  for 
and  j8^: 


ft=  :5-4(rft.)-J^(i-ft)smC 


eBu  dt 


ft 


■;|-4(rft)-;§^(l-ft)cosC. 


(3.34) 

(3.35) 


Substituting  these  equations  into  Eq.  (3.19)  and  making  use  of  the  relations  in  Eqs.  (3.30- 
3.33)  results  in 

A  =  ■  ‘  c  -  ft  c). 

(3.36) 

Using  Eq.  (3.19)  to  solve  for  (j8y  sin  ^  —  jS^cos  ^)  in  terms  of  dy/dt,  Eq.  (3.36)  becomes 


dy  Eo  d  .  ^  Eo  d  .  ^  m  ,  dy 


dt  cB„  dt 


cB„  dt 


eB,, 


dt 


(3.37) 


Extracting  an  overall  time  derivative  gives 


d 

dt 


Eo  „  ^  Eo  „  .  ^  m 

7+  -^rl^ycos  C  +  —yl5xSmC-  -^£COry 

Ct>ii  CiJu  ^t>ii 


=  0. 


(3.38) 


Eactoring  y  and  making  use  of  the  constant  e,  defined  in  Eq.  (3.27),  yields 

=  0 


d 

dt 


J7  77  tyi 

r(  1  +  cos  C  +  -^/3x sin  C  -  — ro(  1  -  ^zo) 

ClJjj  CDii 


(3.39) 
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where  Jo  is  eleetron’s  initial  Lorentz  faetor  value  at  t  =  0,  and  e’s  value  ean  be  set  from 
initial  eonditions  sueh  that  e  =  7^(1  —  Multiplying  both  sides  by  eBu/myo,  we  ean 
integrate  exaetly, 


r 


pJ7 

^  1  -  sin  C  +  cos  C) 

Jotn  Yotnc 


=  eonstant 


In  Eq.  (3.40)  we  ean  define  the  following  eonstants, 

eBu 


(Or 


Join 

eEo 


(Oe  = 

Yotnc 

A(0  =  (Or  -  (Or{l  -  Pzo) 


(3.40) 


(3.41) 

(3.42) 

(3.43) 


where  cOc  is  the  relativistie  eyelotron  frequeney,  Aft)  represents  the  differenee  between  the 
relativistie  eyelotron  frequeney  and  the  Doppler  shifted  light  frequeney  [1],  and  (Oe  is  elee- 
trie  field  analogue  to  the  relativistie  eyelotron  frequeney.  With  these  substitutions,  Eq. 
(3.40)  beeomes 

7[Aft)  +  ft)£(/3xSin^+j8y0os^)]  =  eonstant  =  r/ .  (3.44) 


Two  other  eonstants  of  the  motion  ean  be  developed  from  the  transverse  equations  of  mo¬ 
tion,  Eqs.  (3.20)  and  (3.21).  The  same  proeedure  of  eliminating  (1  —  j5z)cos^  in  favor  of 
a  time  derivative  ean  be  applied  on  Eqs.  (3.20)  and  (3.21).  Eor  example,  substituting  Eq. 
(3.30)  into  (3.20)  yields. 


d 

dt 


(r/3x) 


eEp 

me 


eBu 

m 


Py 


(3.45) 


We  ean  faetor  out  a  time  derivative  and  integrate  the  expression  exaetly  giving  us 


rPx- 


eEo  .  „  ,  eBu  ^  ^ 

- sm  Q  H - y  =  eonstant  =  a. 

me  (Or  me 


(3.46) 
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Similarly,  using  Eq.  (3.31)  in  (3.21),  we  arrive  at 


p]7 

7j8y - ^cos^ - =  constant  =  Z?.  (3.47) 

mcCOr  me 

Table  3.1  summarizes  the  equations  of  motion  and  constants  of  motion  where  the  electron 
phase  is  ^  =  kyZ—  COyt.  Our  goal  is  to  solve  these  equations  for  the  electron’s  position  and 


Table  3.1:  Equations  and  Constants  of  the  Motion 


Eq.  No. 

Differential  Equation 

Eq.  No. 

Constant  of  the  Motion 

3.19 

S  =  -^[^xCOsC-^ySinC] 

3.27 

e  =  r(i-^z) 

3.20 

iir^)  =  -ScosC(i-^.)-^^, 

3.44 

7]  =  Y [^(0 +  (Oe{I3x  sin  ^+I5y  cos  Q] 

3.21 

im)=  ^sinC(l-/3,)  +  ^/3, 

3.46 

3.22 

i(rft)  =  -^i;5,cosc-;SvsmC] 

3.47 

^  =  r^>-i^cosC-^^ 

energy  and  to  determine  the  optimal  conditions  for  significant  energy  gain. 


3.3  Integral  Representation 

In  order  to  solve  Eqs.  (3.19)  through  (3.22)  and  arrive  at  an  expression  that  relates  distance 
traveled  in  the  interaction  region  to  energy  gained,  it  will  be  beneficial  to  recast  some  of 
the  equations  in  Table  3.1.  Consider, 


J3_L-Er  =  (j8y,j8y,0)  ■Eo(cosC,-sinC,0) 
=  E  o[px  cos  i^-py  sin  Q 
=  /3^£’o  cos  9 , 
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where  /3^  =  y  +  /3 ^  and  9{t)  is  the  angle  between  the  electric  field  vector,  Er  and  j3_L 

such  that  0  <  d{t)  <  7t.  Therefore, 


[/3xCos^  —  jSjSin^]  =  j8^cos0.  (3.48) 

Also  consider  the  following  cross-product, 

J3_lX  Er  =  {Px,Py,0)  X  Eo(cosC,-sinC,0) 

=  -Eo  [px  sin  C  +  cos  Qt  (3.49) 

Squaring  both  sides  we  have  and  using  the  magnitude  definition  of  the  cross-product  we 
have, 

(J3^x  Er)-(i5^x  Er)=E^,[l5xCosC  +  PySmC]^  =  PlE^o^m^d.  (3.50) 

Therefore, 

[Px  sin  ^  -h  Py  cos  =  P±  sin  0.  (3.5 1) 


With  Eqs.  (3.48)  and  (3.51)  in  mind,  let’s  return  to  Eqs.  (3.19)  and  (3.44).  We  now  have 


dy 

dt 

V 


~ - P±cosd 

me 

=  7(A(0-f  (O£j8^sin0). 


(3.52) 

(3.53) 


Note  that  dt  =  ydzjciy—  e)  from  Eq.  (3.27).  Substituting  Eqs.  (3.29)  and  (3.42)  into  Eq. 
(3.52),  we  now  have. 


dy  = 


ypOE  [2ye-e^ 
c  (r-e) 


cos  0  dz- 


(3.54) 


We  can  eliminate  the  cos  6  factor  by  solving  Eq.  (3.53)  for  sin0  and  noting  that  sin^  0  + 
cos^  0  =  1,  we  then  have. 


COS0 


± 


to|(27e-e^- 

m|(27e 


l)-(7A(0-r7) 

-e^-  1) 


211/2 


(3.55) 
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Substituting  the  above  equation  in  Eq.  (3.54)  we  arrive  at 


dy=  ±yo 


1 


[cdlilye  -  -  1)  -  (yAco  -  rj)^] 


c(r-e) 


(3.56) 


In  utilizing  the  eonstants  of  the  motion  e  and  rj ,  we  have  arrived  at  an  equation  that  does 
not  explieitly  depend  on  the  angle  6  between  and  Er.  Integrating  both  sides  of  Eq. 
(3.56)  gives  us 


7  -  e 


(o|(27'e  -£2-  1)  _  (7'A(o-t])2 


1/2 


dy' 


(3.57) 


where  the  negative  solution  was  dropped  sinee  we  are  interested  in  energy  gain.  The  left- 
hand  side  is  a  family  of  eurves  known  as  eurtate  eyeloids  [1],  while  the  right-hand  side 
yields  the  distanee  traveled  along  interaetion  region. 


3.4  Curtate  Cycloids 

We  ean  motivate  a  geometrie  understanding  of  Eq.  (3.57)  by  eonsidering  the  mathematies 
of  eyeloids.  A  eyeloid  is  the  eurve  formed  by  traeing  a  point  either  inside,  outside,  or 
on  the  rim  of  a  eirele  of  radius  ^  rolling  along  a  straight  line  [26].  The  general  equation 
eharaeterizing  the  eurve  is  written  in  parametrie  form: 

A  =  (^  t//-  Ksin  t// 

Y  =  ^  —  Kcosxj/  (3.58) 

where  ^  is  the  radius  of  the  eirele,  k  is  the  distanee  from  the  eenter  of  the  eirele  from  where 
the  eurve  is  traeed,  t//  is  the  parametrie  angle.  There  are  three  types  of  eyloids,  as  shown  in 
Eig.  3.2.  Standard  eyeloids  have  k  =  curtate  cycloids  have  k  <  ^,  and  prolate  cycloids 
have  K> 

Consider  a  curtate  cycloid  with  K  <  (^ .  If  we  take  differentials  of  Eq.  (3.58)  then 

dX  -  Kcos  V/)  Jt// 

dY  =K‘sint//Jt//  (3.59) 
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standard  Cycloid 


Standard  Cycloid 


X 


Figure  3.2:  Diagram  describing  different  types  of  cycloids:  standard  cycloid,  curtate,  and  prolate. 

The  radius  of  the  wheel  is  ^ .  The  distance  along  the  radius  where  curve  is  drawn  is 
given  by  fc,  from  [26]. 


Now  consider  the  produet  Y  dY, 

Y dY  =  —  K cos  sin Xjfdxif)  =  KsinXifdX.  (3.60) 


Therefore, 

-  YdY 

dX  =  — ^ - .  (3.61) 

Ksin  t// 

We  can  eliminate  the  sin  y/  term  in  Eq.  (3.61)  with  the  standard  Pythagorean  trigonometric 
identity  and  Eq.  3.58  resulting  in 


dX 


YdY 

^-F^  +  2^F+(k2-<^2) 


(3.62) 
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Notice  this  expression  matehes  the  form  of  Eq.  (3.57).  Integrating  both  sides  of  Eq.  (3.62) 
in  order  to  get  an  explieit  expression  for  x  as  a  function  of  y,  we  have 


K 


(3.63) 


We  have  to  be  eareful  when  employing  the  limits  of  integration  sinee  the  square-root  term 
must  be  real  and  aresine  funetion  is  valid  only  when  the  argument  is  less  than  or  equal  to 
1 1 1 .  The  aresine  requirement  provides  the  following  constraint  on  Eq.  (3.63): 


tzl 

K 


<  1. 


(3.64) 


Squaring  both  sides  of  Eq.  (3.64),  we  have 


-Y^  +  2^Y  +  k^-^^>0. 


(3.65) 


The  radicand  of  the  square-root  term  in  Eq.  (3.63)  matehes  Eq.  (3.65).  Therefore,  only  one 
eonstraint  on  the  value  of  Y  exists.  Solving  the  inequality  gives  the  following  eonstraint  on 

F: 

^-k<Y<^  +  k.  (3.66) 

Eet  us  eonsider  a  simple  example  where  <^  =  3  and  k  =  1  in  Eq.  (3.58).  Therefore,  the 
minimum  and  maximum  values  for  7  are  2  and  4,  respeetively  per  Eq.  (3.66).  The  explicit 
function  for  X  using  Eq.  (3.63)  is 

F(F)  =  -  \/-F^-t-6F-8-3sin^ 

for  2  <Y  <4.  In  order  to  find  T(X)  we  would  need  to 


1(3-7)  +  ^  (3.67) 

invert  Eq.  (3.67)  numerieally  and 
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use  periodic  boundary  conditions  such  that  minima  and  maxima  of  Y  occur  at 

■^max  ~  (2m  +  1)^<3 
Y  max  =  ^  K 

2fmin  ~  2m7r<^ 

Ymin  =  ^-K  (3.68) 


form  =  0, 1,2,3...  [26]. 

In  order  to  make  the  appropriate  comparison  to  Eq  (3.62)  we  must  manipulate  Eq.  (3.57)  a 
bit.  Eet  y  =  y—  e  and  factor  Am  from  the  radicand.  Therefore,  Eq.  (3.57)  becomes 


y'dy 


+  2 


co|e+277Aco-2Aco2e 


y+ 


772+277  Acoe—e^Aa)^ 
A(m2 


Am— z. 


(3.69) 


We  can  now  make  the  following  identifications.  The  y  axis  is  proportional  to  the  energy  of 
electron  and  the  abscissa  is  proportional  to  the  z  distance  traveled  by  the  electron.  Relating 
this  to  the  curtate  cycloid,  the  “radius”  of  the  wheel  is  related  to  the  amplitude  of  energy 
oscillations  as  the  electron  moves  down  the  laser  accelerator.  The  radius  is  given  by 

^  _  m|e  +  2r7Am  -  2^(o^e 

9  ^^2 


and  the  point  from  the  center  of  the  wheel  where  the  curtate  cycloid  curve  is  drawn  is  given 
by 


K  = 


m|e2  _  (ol  _  .^2  IrjAcoe  -  e^Aco^ 


Am2 


+ 


(0‘^e^  +  Ari^Ao)^  +  4Am^e^  +  4m|er]  Am 


4m^Am^e^  -  SriAco^e 


Am4 


(3.71) 


An  explicit  equation  for  distance  as  a  function  of  energy  gained  can  be  immediately  identi¬ 
fied  through  Eq.  (3.63)  with  appropriate  substitutions  for  ^  and  K  given  by  Eqs.  (3.70)  and 
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(3.71)  with  the  constraints  on  7  such  that  ^  —  k<y— e<^  +  K  with  the  form 


Z  oc 


—  7'2  +  2(^  7'  +  2  _  ^ 


K 


(3.72) 


3.5  Resonance  Condition 


The  resonance  condition  is  defined  as 


Am  =  0 

Am=  — -mTl-/3,o)  =0  (3.73) 

7om 


This  is  the  case  where  the  electron’s  relativistic  cyclotron  frequency  is  equal  to  the  Doppler 
shifted  light  frequency.  When  these  two  frequencies  are  balanced,  the  electron  experiences 
a  sustained  acceleration.  Depending  on  initial  conditions,  the  electron  may  initially  gain  or 
lose  energy.  Furthermore,  if  resonance  condition  satisfied,  then  it  is  maintained  throughout 
the  interaction.  Since  £  is  a  constant  of  the  motion,  we  can  rewrite  Eq.  (3.73)  as 


eBu  £ 

- (Or-  =  0. 

7m  7 


(3.74) 


Multiplying  by  yon  both  sides  of  Eq.  (3.74),  we  have 


eBu 

- m,£  =  0.  (3.75) 

m 


Everything  on  the  left-hand- side  of  Eq.  (3.75)  is  a  constant,  therefore  Am  is  zero  for  all 
time. 


There  are  infinitely  many  combinations  of  Jo,  B^,  Xy,  and  /3^o  in  Eq.  (3.43)  that  would  drive 
Am  to  be  zero.  We  can  use  the  resonance  condition  as  a  constraint  to  help  us  sample  the 
parameter  space. 
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3.6  Electron  Phase 


Energy  gained  or  lost  by  electrons  in  the  laser  accelerator  is  given  by  Eq.  (3.52).  Recall 
that  0,  a  dynamical  variable,  is  the  angle  between  j3j_  and  Er  as  shown  in  Eig.  3.1.  If 
6  =  71,  then  the  right-hand  side  of  Eq.  (3.52)  is  positive  since  cos(;r)  =  —  1.  Therefore,  7 
will  grow  with  time.  If  0  =  ;r/2,  the  right-hand  side  is  zero,  7  will  be  at  an  extremum.  The 
seemingly  worst  case  for  a  laser  accelerator  system  is  when  0=0,  then  the  right-hand  side 
of  Eq.  (3.52)  will  have  the  greatest  negative  value  resulting  in  electrons  losing  energy.  The 
initial  angle  between  j3  and  Er  at  the  start  of  the  interaction  is  defined  to  be  60  ■ 

The  phase  for  each  electron  given  by  ^  =  krZ  —  (Ort  will  have  a  value  an  initial  value  of  = 
krZo  at  the  start  of  the  interaction  since  t  =  0.  Experimentally,  the  laser  accelerator  would 
be  accepting  relativistic  electrons  from  an  injector  system.  There  will  be  a  distribution  of 
electrons  spread  along  z  as  shown  in  Eig.  3.3.  This  spread  in  is  equivalent  to  0  <  0o  <  2;r 
over  a  wavelength  of  light.  So  we  expect  some  electrons  to  gain  energy;  some  electrons 
to  lose  energy.  But  as  we  will  see  in  Chapter  5,  having  an  unfavorable  60  may  not  be  so 
detrimental  for  achieving  good  electron  energy  gain  over  the  interaction  length. 


Figure  3.3:  Diagram  of  initial  electron  distribution  in  the  longitudinal  direction  over  two  wave¬ 
lengths  of  light.  In  reality  there  would  be  millions  of  electrons  in  the  distribution,  af¬ 
ter  [27]. 
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3.7  Maximum  Energy  Gain 


The  ideal  case  for  maximum  electron  acceleration  and  energy  gain  is  when  the  resonance 
condition  is  met  or  when  Am  =  0  and  with  Oo  =  Tt.  With  this  substitution,  Eq.  (3.57) 
becomes, 


7 


dy'  =  ^. 


(3.76) 


Jo 


c 


Integrating  both  sides  and  solving  for  z  we  have. 


z  = 


To  3m|e2  [ 

For  large  values  of  y,z°^  y^^^  [1]. 


lai^ey'  -  co^e^  -  m|  -  |  ©ley'  +  (o|e^  +  (o|  +  -  3m|e^ 


J  Jo 

{3.11) 


As  mentioned  before,  there  are  infinitely  many  combinations  of  7o,  B^,  and  for 
resonance  to  be  met.  Practically,  an  experiment  can  be  setup  by  choosing  the  laser  and  an 
external  magnetic  field  leaving  y,  and  as  free  variables  in  order  to  find  resonance.  At 
resonance  using  Eq.  (3.41)  we  have. 


Am 


eBu 

myo 


^zo) 


0 


so, 

pB 

^-^  =  co,{l~M 

myo 


Using  Eq.  (3.27)  we  have 


eBu 

(Ortn 


(3.78) 


With  Eqs.  (3.24)  and  (3.78),  noting  that  =  (y,  —  e)/yo  again  from  Eq.  (3.27),  we  can 
solve  for  the  injection  angle 


a  =  cos 


yo-£ 


=  cos 


To 


eBu 

(Orin 


(3.79) 


In  order  for  a  to  be  a  real-valued  angle,  there’s  a  requirement  on  the  inverse-cosine  function 


39 


of  Eq.  (3.79)  such  that  the  argument  must  be  less  than  or  equal  to  one, 


Yo-e 

VYo-'^ 


<  1. 


(3.80) 


Squaring  both  sides  and  completing  the  inequality  for  fo  gives  us  a  minimum  to  achieve 


resonance, 


7min 


e^  +  1 
2e 


(3.81) 
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CHAPTER  4: 
Numerical  Methodology 


4.1  Finite-difference  Method 


Both  Eqs.  (3.72)  and  (3.77)  represent  analytieal  solutions  for  the  off-resonanee  and  reso- 
nanee  eases,  of  distance  along  z  as  a  function  of  7  respectively;  but  they  are  quite  unwieldy. 
Finite  difference  numerical  integration  methods  provide  a  much  cleaner  solution. 

We  can  rewrite  Eqs.  (3.19)  through  (3.22)  in  terms  of  four-momentum  components.  Re¬ 
calling  that  =  m'Y{c,Vx,Vy,v^  we  have, 


dp^ 

dt 

dp^ 

dt 

dp^ 

dt 

djp 

dt 


[p^cos^  — sin 


-eEo  1 


^  eBu  2 

COSL - P 

ym 


eEo\  1-^  sinC  +  — p 
\  ymc  I  ym 

Tp^cos^  — p^sin^l 

ymc 


(4.1) 

(4.2) 

(4.3) 

(4.4) 


where  ^  =  krZ  —  cod,  =  ymc,  p^  =  yrnfi^c,  p^  =  ym^yC,  and  p^  =  ym^^c.  Turning  the 
first  equation  into  a  finite-difference  equation,  we  have 


7’oid  + 


- [Px  cos  C  -  Pv  sin  At 

ymc  J 


The  “old”  and  “new”  labels  corresponds  to  previous  time  step  values  and  updated  time  step 
values,  respectively,  of  the  variable  in  question.  We  can  then  determine  the  value  of  7new 


through 


7new 


Pnew 

me 


(4.6) 
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Then  using  Eq.  (3.27) 


i8 


new 


Tnew  £ 
Tnew 


since  £  =  7(1  —  is  a  constant  of  the  motion.  The  position  in  z  is  then 


(4.7) 


Znew  ~  Zold  + 


(4.8) 


Similarly  for  the  x  component, 


P 


1 

new 


7’old  + 


-eEo 


cos  ^ 


eBu  9I 
- p  At 

\  old 


Then 


Therefore  the  x-position  is 


^new 


Pnew 

Tiew^C 


^new  ~  ^old  4“  Px  , 


(4.9) 


(4.10) 

(4.11) 


and  similarly  for  the  y  component.  Since  we  are  updating  the  new  positions  on  the  updated 
velocities,  this  is  known  as  Euler-Cromer  method  [28]. 


Numerical  integration  using  a  4th  order  Runga-Kutta  (RK4)  method  was  also  performed. 
The  percent  difference  in  electron  energy  gain  between  RK4  and  Euler-Cromer  methods  is 
only  0.03  percent.  Since  the  percent  difference  is  negligible,  Euler-Cromer  methods  were 
employed  in  the  dissertation  due  to  their  speed  and  simplicity. 


4.2  Initial  Conditions  on  Resonance  (Ao)  =  0) 

The  task  in  this  section  is  to  numerically  solve  the  equations  in  Table  3.1  subject  to  the 
electrons  entering  the  interaction  region  with  the  resonance  condition  satisfied  such  that 
Am  =  0.  We  ultimately  want  to  know  the  energy  gained  by  the  electron  and  its  trajectory. 

Assume  that  the  initial  electron  energy  and  external  magnetic  field  are  independently  cho¬ 
sen  parameters.  Erom  Eq.  (3.79),  we  have  an  equation  for  the  injection  angle,  a  = 
a(7o,fi„),  as  a  function  of  two  free  variables.  The  following  result  shown  in  Eig.  4.1  is 
a  surface  of  required  injection  angles  for  resonance. 
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If  we  select  a  magnetic  field  strength  of  45  T,  then  we  can  slice  through  a  constant  external 
magnetic  field  in  Fig.  4.1  resulting  in  Fig.  4.2.  We  now  have  an  optimal  injection  angle 
for  various  initial  electron  energy  values  where  the  minimum  value  of  Jo  in  Fig.  (4.2)  is 
consistent  with  Eq.  (3.81). 
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Required  injection  angle  for  =  45  T 


Figure  4.2:  Shown  is  the  allowed  ordered  pair  combination  of  initial  electron  energy,  and  initial 
injection  angle,  a  for  a  given  magnetic  field  to  guarantee  resonance  for  a  given  (Or- 


Consider  the  eleetron  at  the  start  of  the  laser  aeeelerator  system  at  time,  t  =  0.  Using  Eqs. 
(3.15)  and  (3.17),  the  eleetrie  field,  Er,  is  then  oriented  along  the  jc  axis  or  Er  =  EgX  sinee 
=  0)  =  krZ  —  COrt  =  0-  From  Eq.  (3.52)  it  would  seem  that  the  greatest  energy  gain 
would  require  an  anti- alignment  between  Er  and  sueh  that  6  =  n  initially.  The  other 

spherieal  eoordinate  angle,  ^o,  whieh  is  the  angle  between  and  is  then  also  equal 
to  n.  In  general,  may  take  on  any  value  between  0  and  n.  So,  with  a  and  established, 
the  initial  eonditions  for  position  and  veloeity  ean  be  ealeulated  using  Eqs.  (3.23)  through 
(3.26). 

4.3  Electron  Trajectory 

For  plane- wave  propagation  (no  diffraetion)  with  6  =  n  for  maximum  aeeeleration,  both 
Eqs.  (3.52)  and  (3.56)  are  equivalent  with  Aco  =  T]  =  0.  It  is  straightforward  to  numerieally 
integrate  and  update  7.  Then  for  eaeh  timestep  using  the  eonstant  of  the  motion  e,  /3^  ean 
be  updated  and  the  position  z  ealeulated.  The  values  for  Px  and  j8y  ean  also  be  evaluated  by 
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making  use  of  constants  a  and  b  where  their  values  are  given  by  initial  conditions.  From, 
Eqs.  (3.46)  and  (3.47),  and  jiy  are 


/3x 

Py 


a  H - sm  Q - y 

mcCOr  me 

b  H - cos  g  H - .r 

me  (Or  me 


(4.12) 

(4.13) 


We  can  determine  the  number  of  times  an  electron  will  spiral  around  the  external  magnet’s 
field  lines  by  integrating  Eq.  (3.41)  (and  making  the  substitution  of  y{t)  for  y^)  over  the 
time  it  takes  the  electron  to  clear  the  interaction  region.  Stronger  magnetic  fields  results 
in  tighter  and  more  rotations.  Trajectories  from  numerical  integration  will  be  shown  in 
Chapter  5. 


4.4  Synchrotron  Radiation 


An  estimate  of  power  loss  through  synchrotron  radiation  can  be  made  if  we  consider  the 
simplified  case  with  the  laser  source  switched  off.  With  only  the  external  magnetic  field 
and  no  laser  field  present,  Eo  =  0,  Eqs.  (3.19)  through  (3.22)  reduces  to 


y{t)  =  7(0)  =  7o,  constant 

dvx  _  eBu 
dt  my  ^ 

dvy  eBu 

EJE  =  - 

dt  my 


=  ^z(O)  =  Vzo,  constant 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


Eqs.  (4.15)  and  (4.16)  can  be  solved  by  letting  p  =  Vx  +  ivy  and  (Oc  =  eBu! my.  Therefore, 
p  =  icOep.  Solving  this  with  initial  conditions  yields  p(t)  =  where  v^o  = 


45 


vlg  +  and  5  =  tan  ^  {vyo/vxo)  ■  Integrating  once  more  yields  the  following 


x{t)  = 
y{t)  = 
zit)  = 


■^^sinfcOct  +  5)  +Xo  —  -^^sinfS) 

COc  ®c 

-■^^cos(tt)ct  +  5)  +yo  +  ^^cos(5) 

(Oc 

Zo  “h  . 


(4.18) 


The  radius  of  the  orbit  is  then  given  by  the  cyclotron  radius  for  an  electron  in  a  static 
magnetic  field, 


r=  \j{x{t)-Xof  +  {y{t)-yof 
^  ymv^o 
eBu 


(4.19) 


We  can  calculate  the  synchrotron  power  loss  using  the  relativistic  Larmor  formula  given 
by  [29], 


P  = 


Boe^fc 

6it 


/S"-  jSxjS 


(4.20) 


Taking  one  and  two  time  derivatives  of  Eqs.  (4.18)  and  substituting  into  Eq.  (4.20)  yields. 


67t 


P  = 


6;r 


We  arrive  at  [29], 


p_Poffc  202 


(4.21) 


Plugging  in  nominal  values  for  an  external  magnetic  field  of  45  T  and  300  T,  we  can  expect 
synchrotron  energy  losses  to  be  on  the  order  of  several  electron-volts  (eV)  to  5  keV,  respec¬ 
tively.  These  figures  represent  an  upper-bound  estimate  energy  losses.  Once  the  laser  field 
is  present,  as  will  be  shown  later,  7  increases  however  /3^  is  small  and  its  value  decreases. 
Eor  energy  gain  of  approximately  one  GeV,  power  loss  through  synchrotron  radiation  is 
negligible. 
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4.5  Numerical  Integration 

The  numerical  integration  procedure  for  off-resonance  or  resonant  analysis  begins  by  se¬ 
lecting  global  constants  such  as  laser  wavelength,  laser  power,  estimated  spot  size,  external 
magnetic  field  strength,  and  time  steps.  If  resonance  analysis  is  required,  the  angle  for  res¬ 
onance  based  on  the  global  constants  is  calculated.  We  compute  initial  momentum  values 
in  preparation  for  numerical  integration.  Based  on  the  initial  values,  we  perform  Euler- 
Cromer  integration  to  determine  “new”  values  for  7  and  the  momenta.  We  use  these  new 
values  to  calculate  position  and  the  overall  trajectory  of  the  electron.  We  also  calculate 
power  loss  through  synchrotron  radiation  and  perform  diagnostic  checks  of  the  code  for 
each  time  step.  Finally,  we  keep  track  of  the  values  by  writing  the  data  to  a  file  for  post¬ 
processing  in  MatLab. 


47 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


48 


CHAPTER  5: 

Numerical  Analysis  Results 


In  this  chapter,  we  will  use  the  numerical  integration  method  previously  described  to  obtain 
simulation  results  on  the  laser  accelerator  performance  for  a  variety  of  initial  conditions 
including:  on  resonance,  different  values  of  the  initial  phase  0^,  the  angle  between  J3  and 
Er,  off  resonance,  and  real  field  effects. 


5.1  On-resonance  with  6o  =  7t 

Recall  that  resonance  is  defined  as  Am  =  0  and  occurs  when  the  electron  “feels”  the  Doppler 
shifted  circularly-polarized  laser  light  oscillating  at  the  same  frequency  as  the  electron’s 
spiraling  motion  in  the  uniform  magnetic  field.  In  the  simulation,  we  assume  that  the 
external  magnetic  field  and  laser  source  parameters  (laser  power  and  wavelength)  are 
fixed.  The  resonance  condition  is  then  a  function  of  the  injection  angle  a  and  injection 
energy  as  shown  in  Fig.  4.2. 

Let  us  assume  the  following  for  the  simulation  of  a  single  electron  .  The  laser  source  will 
be  like  the  Texas  Petawatt  Laser  (TPL)  system.  The  laser  Rayleigh  length,  Z^,  is  assumed 
to  be  longer  than  interaction  region  of  one  meter  in  order  to  eliminate  diffraction  effects. 
Gaussian  beam  analysis  will  be  considered  later  in  the  chapter.  The  optical  mode  waist  is 
about  1  mm  throughout  the  interaction  region.  The  external  magnetic  field  is  assumed  to 
be  static,  uniform  field  throughout  the  interaction  region.  Strong  magnetic  fields  support 
electron  to  light  coupling  for  energy  gain.  Let  us  also  assume  6o  =  7t  so  that  is  anti¬ 
parallel  to  the  electric  field  of  the  light  as  shown  in  Fig.  3.1. 

For  the  single  electron  simulation.  Table  5.1  summarizes  the  input  parameters. 

Performing  the  numerical  integration  as  described  in  the  previous  chapter,  we  obtain  the 
results  shown  in  Fig.  5.1.  The  results  are  compared  to  the  analytic  solution  given  in  Eq. 
(3.77).  The  percent  difference  is  about  0.2  percent  for  the  final  energy  of  the  electron  after 
a  one  meter  long  interaction.  We  compared  multiple  simulation  results  to  the  analytic  pre¬ 
diction  in  regimes  where  the  analytic  formation  remained  valid.  In  all  cases,  the  simulation 
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Table  5.1:  Input  Parameters  for  Resonance  Laser  Accelerator  Simulation 


Variable 

Symbol 

Value 

Initial  Electron  Lorentz  Factor 

7o 

119 

Initial  Electron  Energy 

So 

61  MeV 

Easer  Peak  Power 

P 

1  PW 

Easer  Wavelength 

Xr 

1  /im 

External  Magnetic  Field 

45  T 

First  Constant  of  Motion 

£ 

4.2  X  10-3 

Initial  Azimuthal  Angle 

Bo 

n  radians 

Injection  Angle 

a 

20.8  /i  radians 

Time  Step 

At 

50  fs 

always  closely  matched  the  analytic  prediction. 


z  (m) 


Figure  5.1:  Comparison  of  analytic  solution  to  numerical  simulation  results  of  an  electron  in  the 
laser  accelerator.  Energy  gain  versus  distance  is  plotted.  The  black  curve  describes  the 
analytic  solution  and  the  red  curve  describes  the  numerical  solution.  The  difference  is 
about  0.2  percent  for  the  final  energy  gain,  which  is  approximately  1.17  GeV. 
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Now  that  we  have  established  good  agreement  between  analytic  and  numerical  solutions 
we  can  now  explore  how  other  parameters  such  as  electron  radial  position  and  velocity 
evolution  over  the  interaction  region.  Fig.  5.2  shows  simulation  results  for  two  electrons 
at  two  different  external  magnetic  field  values,  45  T  (blue)  and  100  T  (red)  with  injection 
angles  a  of  20.8  /irad  and  9.3  mrad  respectively  for  resonance  as  per  Eq.  (3.79).  The  other 
parameters  are  the  same  as  in  Table  ??. 


z  (m)  z  (m) 


Figure  5.2:  Electron  simulation  results  for  the  parameters  given  in  Table  5.1.  The  blue  curve  give 
results  for  B„  =  45  T  while  the  red  curve  give  results  for  B„  =  100  T.  The  top  left 
graph  plots  electron  energy  versus  distance  traveled.  The  top  right  graph  plots  the  radial 
distance  of  the  electron  from  the  optical  axis.  The  bottom  left  graph  plots  [5±  versus 
distance  traveled.  The  bottom  right  graph  plots  1  —  versus  distance  traveled. 


The  top  left  curve  shows  that  the  electron’s  energy  increases  to  nearly  1 .2  GeV  for  45  T 
and  1.6  GeV  for  100  T  in  one  meter.  The  electron’s  orbit  spirals  out  to  a  radius  of  about 
0.17  mm  for  45  T  and  0.15  mm  for  100  T  as  shown  in  the  top  right  graph.  The  value  of 
1  —  /3^  is  plotted  in  the  bottom  right  curve  since  {5^  is  nearly  unity  for  the  entire  length  of  the 
interaction  region.  Note  that  (5^  increases  while  decreases  during  most  of  the  interaction. 
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In  the  bottom  left  curve,  we  can  see  that  j5±  has  a  peak  then  slowly  decreases  for  the  45  T 
case  while  the  100  T  does  not  have  a  peak.  Let  us  try  to  understand  why  has  a  peak  for 
the  45  T  case  while  no  peak  exists  for  the  100  T  case.  It  is  not  immediately  clear  from  Eq. 
(3.29)  alone  that  /3^  should  have  a  maximum.  Calculating  its  derivative  we  have 


r  1 

dz  dz  [fiz) 

-Sa 


27(z)e-e2-  1 


[2riz)£-£^-l] 


2  11-1/2  1  .2  ill/2l^/r(z) 


Setting  Eq.  (5.1)  equal  to  zero  we  have  an  extremum  when 

*  e^+ 1 

r  =  — 


(5.1) 


(5.2) 


Notice  that  y*  =  lymin  from  Eq.  3.81.  Calculating  y*  for  45  T  and  100  T  cases  we  have 


7^(45  T)  ^238 
7* (lOOT)  ^  107 


The  value  for  7*  (100  T)  explains  why  there  is  no  peak  for  the  100  T  case  in  Eig.  5.2  since 
the  simulation  began  with  an  initial  7^  =  119  T.  Substituting  Eq.  (5.2)  back  into  Eq.  (3.29) 
yields 

^±,max  =  ,  ^  =;  (5.3) 

Ve^Tl 

which  agrees  with  the  simulation  value  for  j8^,max-  In  Eig.  5.3,  we  have  again  plotted  the 
electron  trajectories  for  two  the  cases:  45  T  field  in  the  blue  curve  and  100  T  field  in  the 
red  curve.  Notice  as  7  increases  along  the  z  axis  it  causes  the  radius  of  spiraling  motion 
around  the  magnetic  field  line  to  increase.  There  are  more  numerous  spirals  for  the  stronger 
field.  The  number  of  spirals  gives  an  approximate  number  of  wavelengths  of  light  that  pass 
over  the  electron  near  resonance.  High  powered  lasers  such  as  the  Texas  Petawatt  Easer 
produces  Gaussian  pulses.  The  Eourier  transform  of  these  pulses  would  also  have  a  spread 
of  frequencies  in  the  frequency  domain.  As  shown  in  Eig.  5.3,  only  about  ten  wavelengths 
of  light  pass  over  the  electron.  To  first  order  we  do  not  need  to  be  concerned  with  the  effects 
to  the  frequency  changes  associated  with  the  pulse  shape. 
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Figure  5.3:  Electron  trajectories  in  the  laser  accelerator  for  two  different  external  magnetic  Field 
CasesLaser  accelerator  for  the  parameters  given  in  Table  5.1  in  the  blue  curve.  The  100 
T  field  trajectory  is  shown  in  the  red  curve  for  comparison. 


Table  5.2  summarizes  the  simulation  results  for  both  eases. 

Table  5.2:  Results  for  Single  Electron  Simulation  of  Easer  Accelerator  at  45  T  and  100  T 


Variable 

Symbol 

Value  at  45  T 

Value  at  100  T 

Final  Lorentz  Factor 

Yf 

2282 

3109 

Final  Energy 

1.17  GeV 

1.59  GeV 

Final  Electron  radial  position 

r 

0.16  mm 

0.13  mm 

5.2  Electron  Phase  Performance 

A  realistic  electron  bunch  is  thousands  of  times  longer  than  a  wavelength  of  light.  As 
mentioned  in  §  3.6,  there  will  be  millions  of  electrons  in  each  optical  wavelength  that  are 
randomly  distributed  in  the  microscopic  longitudinal  direction.  Recall  that  ^  =  kyZ  —  (Oft 
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and  at  time  t  =  0,  l^o  =  hzo  =  liiZolK-  Their  initial  ^o’s  will  therefore  span  a  1%  range 
over  a  distance  where  Zo  changes  by  a  wavelength  of  light,  Xy.  We  expect  some  electrons 
to  gain  energy  and  some  electrons  to  lose  energy  depending  on  ^o- 

We  select  the  coordinate  system  so  that  electron  is  initially  at  the  origin  and  the  laser  electric 
field  is  initially  along  the  a:-axis.  In  order  to  sample  different  microscopic  values  of  l^o,  we 
consider  three  representative  values  of  =  0,  ;r/2,  ;r  in  the  simulations.  Recall  from  Eq. 
(3.52),  that  these  values  are  anticipated  to  cause  the  electron  to  initially  lose  or  gain  energy. 

Consider  the  following  simulation,  shown  in  Fig.  5.4  with  =  130  (above  the  minimum 
energy  required  for  resonance),  =  45  T,  and  on  resonance  injection  with  Ores  ~  2.4  mrad. 

In  Fig.  5.4,  we  examine  the  evolution  of  energy  versus  distance  traveled  for  the  first  15mm. 
The  middle  graph  shows  versus  distance.  The  bottom  graph  shows  the  evolution  of 
9 .  The  colored  curves  shows  the  evolution  of  the  respective  parameter  for  three  different 
initial  9o.  The  black  curves  shows  9o  =  K,  the  green  curves  show  9o  =  njl,  red  curves 
show  9o  =  0. 


In  describing  the  energy  plots  we  must  recall  Eq.  (3.52) 


dy  eEo 

—  = - p±cos0. 

at  me 


Notice  the  factor  of  cos  6  on  the  right-hand  side  of  the  energy  equation.  For  the  do  =  it 
case  shown  by  the  black  curve,  the  energy  initially  increases  since  the  right-hand  side  of 
Eq.  (3.52)  is  positive.  For  the  the  =  ;r/2  case  shown  by  the  green  curve,  the  energy  is 
initially  constant  since  the  right-hand  side  of  Eq.  (3.52)  is  zero.  For  the  Bo  =  0  case  shown 
by  the  red  curve,  the  energy  initially  decreases  since  the  right-hand  side  of  Eq.  (3.52)  is 
negative. 


Fet  us  now  look  at  the  perpendicular  velocity  graphs  for  j5±.  For  both  the  Bq  =  n  and 
Bo  =  7r/2  cases,  /3^  goes  up  initially,  as  shown  by  the  black  and  green  curves,  consistent 
with  the  electron’s  increasing  energy  from  the  top  graph.  For  the  Bo  =  0  case  shown  by 
the  red  curve,  /3j^  decreases  in  value  until  it  reaches  zero.  When  I5±  goes  to  zero,  the  angle 
B,  shown  in  the  bottom  graph,  changes  from  0  to  ;r  quickly.  Then  /3^  and  the  energy  both 
begin  to  grow.  This  is  a  surprising  result.  It  was  not  expected  that  an  electron  that  starts  out 
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z  (mm) 

Figure  5.4:  Simulation  results  for  on-resonance  but  different  initial  do  for  jo  =  130,  B„  =  45.  The 
top  graph  shows  energy  versus  distance  traveled  for  the  first  15mm.  The  middle  graph 
shows  versus  distance.  The  bottom  graph  shows  the  evolution  of  6.  The  black  curves 
shows  do  =  71.  The  green  curves  show  do  =  Tl/l.  Red  curves  show  do  =  0. 


moving  parallel  to  the  light’s  eleetrie  field  would  eventually  gain  energy. 

Let  us  consider  a  similar  analysis  of  Of,  but  with  lower  initial  electron  energy  of  Jo  =  125 
and  Ofres  ~  1.8  mrad.  These  results  are  shown  in  Fig.  5.5.  The  electron’s  perpendicular 
velocity  /3^  starts  out  with  a  smaller  value  than  the  Jo  =  130  case  since  the  required  injec¬ 
tion  angle  is  smaller.  So  for  the  0o  =  0  case,  /3^  goes  to  zero  sooner  along  the  z  axis.  This 
causes  6  to  change  to  sooner  than  the  Jo  =  130  case.  We  again  see  that  an  electron  that 
starts  the  interaction  with  an  unfavorable  6o,  eventually  changes  to  a  favorable  0.  Thus 
they  accelerate  to  high  energies  along  with  electrons  that  initially  started  with  favorable 
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Figure  5.5:  Simulation  results  for  on-resonance  but  different  initial  do  for  jo  =  125,  B„  =  45.  The 
top  graph  shows  energy  versus  distance  traveled  for  the  first  15mm.  The  middle  graph 
shows  versus  distance.  The  bottom  graph  shows  the  evolution  of  6.  The  black  curves 
shows  do  =  71.  The  green  curves  show  dg  =  7i/2.  Red  curves  show  do  =  0. 

do  values.  Again,  this  is  a  surprising  and  encouraging  result.  Initially,  it  was  thought  that 
electrons  that  did  not  have  a  favorable  do  would  not  gain  energy. 

Let  us  again  consider  a  similar  analysis  of  initial  do  but  with  an  even  lower  initial  electron 
energy  of  Jo  =  120  and  ctres  ~  Imrad.  These  results  are  shown  in  Fig.  5.6.  With  =  120 
case,  /3^  starts  out  with  an  even  smaller  value  than  the  previous  two  cases  since  the  required 
injection  angle  is  smaller.  Therefore,  /3^  goes  to  zero  sooner  which  causes  do  to  change 
quickly  7t  sooner. 

Let  us  again  consider  a  similar  analysis  of  initial  do  but  an  initial  electron  energy  of  7o  =  1 19 
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Figure  5.6:  Simulation  results  for  on-resonance  but  different  initial  do  for  Yo  =  120,  B„  =  45.  The 
top  graph  shows  energy  versus  distance  traveled  for  the  first  15mm.  The  middle  graph 
shows  versus  distance.  The  bottom  graph  shows  the  evolution  of  6.  The  black  curves 
shows  do  =  71.  The  green  curves  show  do  =  Tzjl.  Red  curves  show  do  =  0. 


and  Ctres  ~  20  /frad.  These  results  are  shown  in  Fig.  5.7.  We  are  now  quite  familiar  with 
the  story.  In  the  middle  graphs,  /3^  starts  with  a  mueh  smaller  value  than  the  previous  three 
initial  energy  cases.  For  the  0o  =  0  case,  /3^  evolves  to  zero  even  sooner  which  quickly 
leads  to  0  =  ;r.  This  result  encourages  us  to  inject  at  the  minimum  energy  possible  where 
the  resonance  condition  will  remain  satisfied,  given  in  Eq.  (3.81),  if  the  desire  is  for  all  9o 
to  be  accelerated  quickly. 

Let  us  once  more  consider  a  similar  analysis  of  different  9o  but  with  the  lowest  an  initial 
electron  energy  that  will  satisfy  the  resonance  condition.  We  have  Jo  =  7min  and  (Zres  =  0 
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Figure  5.7:  Simulation  results  for  on-resonance  but  different  initial  do  for  Yo  =  119,  =  45.  The 

top  graph  shows  energy  versus  distance  traveled  for  the  first  15mm.  The  middle  graph 
shows  versus  distance.  The  bottom  graph  shows  the  evolution  of  6.  The  black  curves 
shows  do  =  71.  The  green  curves  show  do  =  Tl/l.  Red  curves  show  do  =  0. 


rad.  These  results  are  shown  in  Fig.  5.8.  Minimum  energy  injeetion  results  in  /3^  =  0.  To 
prove  this,  substitute  the  analytie  expression  for  minimum  energy  required  for  resonanee, 
Eq.  (3.81),  into  the  expression  for  /3^  in  terms  of  7  and  £,  Eq.  (3.29).  With  /3^  =  0,  the 
veetor  j3j_  has  no  length.  Therefore,  0  is  undefined.  In  fact,  all  electrons  starting  at  the 
minimum  energy,  will  have  the  same  evolution  for  energy  and  /3^  as  shown  in  Fig.  5.8. 
The  laser  accelerator  for  energies  close  to  the  minimum  required  for  resonance  will  accept 
all  values  of  do  and  resulting  in  energy  gain. 
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Figure  5.8:  Simulation  results  for  on-resonance  with  = /min  ~  118.99927,  =  45.  The  top  graph 

shows  energy  versus  distance  traveled  for  the  first  15mm.  The  bottom  graph  shows 
versus  distance. 


Let  us  now  try  to  understand  why  this  surprising  result  makes  physical  sense.  Recall  that 
this  entire  discussion  has  been  on  resonance  or  when  Am  =  0.  Recall  Eq.  (3.53)  where 
rj  =  /(Am  -F  m^jS^sinO).  Since  we  know  that  is  zero  at  some  point  in  the  evolution 
for  the  Oo  =  0  case,  that  pins  the  value  of  rj  to  zero  for  all  time  since  it  is  a  constant  of  the 
motion.  But  later  in  the  evolution,  /3^  and  /  grow.  The  only  way  to  ensure  that  rj  stays  zero 
is  for  sin0  to  be  zero  or  when  9  =  K.  The  next  line  of  reasoning  eliminates  the  possibility 
that  9  evolves  to  0. 

When  is  zero,  then  and  are  also  zero.  Going  back  to  the  original  equations  of 
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motion,  Eqs.  (3.19)  through  (3.22),  we  now  have 


dt 

d  ,  ^  ^ 

cos  ^ 

=  - 

l-i8. 

me 

eEn  ,  ^  , 

sin^ 

j,m)  = 

me 

4(rft)  ^0 

(5.4) 

(5.5) 

(5.6) 

(5.7) 


So  the  terms  that  deal  with  the  magnetie  field  are  eompletely  turned  off.  They  eontribute 
nothing  when  the  electrons  move  parallel  to  the  field  lines.  Rewriting  Eqs.  (5.5)  and  (5.6) 
with  the  components  of  the  electric  field  factored  out  we  have 


=  (  E^) 

yme^ 

(5.8) 

i(rft) 

=  (  Ey) 

yme  ^ 

(5.9) 

The  right-hand  side  of  Eqs.  (5.8)  and  (5.9)  are  proportional  to  the  force  that  the  electron 
feels  when  it  is  moving  in  a  direction  that  is  opposite  to  the  electric  field  or  when  6  =  n. 

In  this  section  we  have  considered  the  performance  of  the  laser  accelerator  for  only  three 
Oo  cases.  But  we  are  encouraged  that  even  with  the  most  unfavorable  initial  6o,  energy  gain 
is  possible  if  the  initial  energies  are  near  the  minimum  required  for  resonance. 


5.3  Off-resonance  Performance 

As  already  discussed,  the  electron  bunch  in  an  actual  experiment  is  much  longer  than  the 
optical  wavelength  and  therefore  spans  over  the  entire  9o  range  from  0  to  2;r.  Additionally, 
electrons  will  not  be  perfectly  collimated  and  will  consequently  have  a  finite  spread  in 
injection  angles  a.  This  latter  effect  implies  that  the  resonance  condition.  Am  =  0,  will  not 
be  satisfied  for  all  electrons  in  the  bunch  since  for  a  given  initial  energy,  magnetic  field, 
and  laser  frequency.  Am  is  a  function  of  the  injection  angle  a.  Substituting  Eqs.  (3.24)  and 
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(3.41)  in  (3.43)  we  have 


^  «)  ■  (5.10) 

The  right-hand  side  of  Eq.  (5.10)  equals  zero  when  a  =  Ofres  (whieh  is  ealeulated  in  Eq. 
(3.79)).  We  want  to  explore  the  behavior  of  the  laser  accelerator  for  off-resonance  cases  for 
a  spread  of  a  near  Ores- 

Eet  us  begin  with  the  case  7o  =  130  with  an  external  field  'Qu  =  45  T.  The  results  are  shown 
in  Eig.  5.9  where  we  have  conducted  many  simulations  of  the  laser  accelerator  over  a  meter 
long  interaction  and  plotted  the  final  energy  values  Jf  as  a  function  of  a  and  do-  Notice 
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Figure  5.9:  Surface  graph  of  final  electron  energy  jf  versus  initial  do  and  a  for  7,  =  130  and  Bj, 
45.  There  is  a  peak  energy  gain  near  oCres  ~  2.4  mrad  and  do  =  tt  as  expected. 


that  there  is  a  peak  energy  gain  near  Ores  ~  2.4  mrad  and  do  =  TT  as  expected. 

We  can  get  a  better  perspective  by  considering  the  contour  plot  shown  in  Eig.  5.10.  The 
results  indicate  that  electrons  with  initial  conditions  inside  the  blue  region  achieve  99  per- 
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Figure  5. 10:  Contour  graph  of  electron  energy  versus  do  and  injection  angles  a  for  7o  =  130  and  B„ 
=  45.  The  peak  occurs  near  Ores  ~  2.4  mrad.  Electrons  with  do  and  a  inside  the  blue 
region  achieve  99  percent  of  maximum  energy  gain.  Electrons  with  do  and  a  inside 
the  green  region  achieve  99. 1  percent  of  maximum  energy  gain.  Electrons  with  do  and 
a  inside  the  maroon  region  achieve  99.9  percent  of  maximum  energy  gain. 


cent  of  maximum  energy  gain.  Eleetrons  with  initial  conditions  inside  the  green  region 
aehieve  99.1  pereent  of  maximum  energy  gain.  Eleetrons  with  initial  eonditions  inside  the 
red  region  aehieve  99.9  pereent  of  maximum  energy  gain.  These  results  give  a  sense  of 
the  energy  spread  of  the  laser  aeeelerator  output  beam,  depending  on  the  alignment  of  the 
input  beam  from  the  eleetron  injeetor  and  the  spread  in  injeetion  angles. 

Eet  us  consider  a  similar  plot  but  now  for  the  case  of  =  120  with  an  external  field  of 
45  T  shown  in  Eig.  5.11.  These  results  look  similar  to  the  previous  eontour  plot  but  the 
parameter  space  area  resulting  in  99  pereent  energy  gain  has  increased  dramatieally  to 
inelude  all  initial  Oo  values.  Erom  the  diseussion  on  phase  from  §  5.2,  this  makes  sense. 
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Figure  5.11:  Contour  graph  of  electron  energy  versus  6o  and  injection  angles  a  for  7o  =  120  and  B„ 
=  45.  The  peak  occurs  near  Ores  ~  1-8  mrad.  Electrons  with  do  and  a  inside  the  blue 
region  achieve  99  percent  of  maximum  energy  gain.  Electrons  with  do  and  a  inside 
the  green  region  achieve  99. 1  percent  of  maximum  energy  gain.  Electrons  with  do  and 
a  inside  the  maroon  region  achieve  99.9  percent  of  maximum  energy  gain. 


Electrons  with  an  unfavorable  do,  lose  since  j5±  begins  with  a  smaller  value.  After 
goes  to  zero,  The  value  of  6  is  no  longer  meaningful  and  can  change  to  n  for  the  reasons 
discussed  previously.  Now  these  electrons  now  are  able  to  achieve  energy  gain. 

Let  us  consider  a  similar  plot  but  now  for  the  case  of  7o  =  1 19  with  an  external  field  of  45 
T  shown  in  Fig.  5.12.  The  parameter  space  area  resulting  in  99.1  percent  energy  gain  has 
increased  dramatically  to  include  all  initial  6o-  Notice  that  for  Figs.  5.10,  5.11,  and  5.12, 
we  gained  greater  acceptance  range  in  6o  at  the  expense  of  the  maximum  energy  jf  since 
it  went  down  from  Yf  =  2056  to  7/  =  2042. 
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Maximum  y  =  2042 


o  (mrad) 

Figure  5.12:  Contour  graph  of  electron  energy  versus  do  and  injection  angles  a  for  7o  =  120  and  B„ 
=  45.  The  peak  occurs  near  Ores  ~  20  /trad.  Electrons  with  do  and  a  inside  the  blue 
region  achieve  99  percent  of  maximum  energy  gain.  Electrons  with  do  and  a  inside 
the  green  region  achieve  99. 1  percent  of  maximum  energy  gain.  Electrons  with  do  and 
a  inside  the  maroon  region  achieve  99.9  percent  of  maximum  energy  gain. 


Encouraged  by  the  previous  contour  plots  and  the  results  of  §  5.2,  it  seems  that  starting 
with  a  smaller  /3^  ensures  more  electrons  receive  similar  accelerations  and  energy  gains. 
So  let  us  just  start  with  a  beam  entering  on-axis  where  oc  =  /3^  =  0.  Now  Act)  is  only  a 
function  of  two  independent  variables,  Yo  and  We  have  for  on-axis  injeetion 

Acooa  =  —  -  cOr  ( \  -  -  l]  .  (5.11) 

myo  V  ) 

By  varying  yo  and  'Qu  we  expect  the  maximum  energy  to  be  in  a  band  that  follows  the 
Act)  =  0  contour  seen  in  Fig.  5.13. 
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Figure  5.13:  Simulation  results  showing  an  intensity  plot  of  final  electron  energies  with  initial  con¬ 
ditions  that  explore  a  range  of  external  magnetic  fields  from  B„  =  30  T  to  60  T  and 
7o  =  80  to  130.  The  dotted  line  indicates  the  Acoqa  =  0  contour. 


In  Fig.  5.13,  we  have  an  intensity  plot  that  depiets  the  final  eleetron  energies  of  eleetrons 
with  varying  initial  jo  and  external  magnetie  field  values  but  with  on-axis  injeetion  where 
a  =  =  0.  Therefore  the  angle  between  and  Er  has  no  meaning.  The  optieal  foree 

on  the  eleetron  from  the  eleetrie  field  aets  on  the  eleetron  sueh  that  6  quiekly  evolves  to 
a  value  of  n  as  diseussed  at  the  end  of  §  5.2.  The  eleetrons  then  proeeed  to  gain  energy. 
The  resonanee  Act)  =  0  level  eurve  is  overlaid  on  the  intensity  plot  and  shown  as  a  blaek- 
dotted  eurve.  Notiee  how  the  eleetrons  with  good  energy  gain  follows  that  eontour  eurve  as 
expeeted.  There  exists  a  width  about  the  eontour  line  Act)  =  0  for  whereby  the  laser  aeeel- 
erator  performs  well  with  Yf  ^  2500.  This  eorresponds  to  eleetron  energy  gain  of  1.2  GeV 
in  one  meter.  Onee  the  right  eombination  of  parameters  for  resonanee  are  identified,  the 
laser  aeeelerator  is  quite  forgiving  in  its  aeoeptanee  margins  in  initial  energy  and  external 
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magnetic  field  strength  for  on-axis  injection. 

We  can  now  make  an  estimate  on  normalized  emittance  values  that  the  laser  accelerator 
would  accept  from  an  electron  injector.  Normalized  emittance  is  given  by 

£n  =  Y<r><  a>  (5.12) 

Let  us  consider  an  initial  electron  energy  of  7=  130,  <  r  >  =0.5  mm  since  the  optical  waist 
is  about  1  mm,  and  <  a  >  =  0.5  mrad.  The  spread  in  injection  angles  a  is  a  conservative 
estimate  based  on  the  simulations  presented.  We  then  have  =  32.5  mm-mrad,  which  is 
quite  large.  Electron  injectors  with  good  emittance  values  are  about  10  mm-mrad. 


5.4  Magnetic  Field  Effects 

The  external  magnetic  field  source  for  the  laser  accelerator  might  be  constructed  of  resistive 
or  superconducting  solenoids  such  as  those  at  NHMFL.  The  magnetic  field  on  the  axis  for 
a  solenoid  can  be  written  as  [25]: 


oc  5„(cos0i -I-COS02)  (5.13) 

where  0i  and  62  are  defined  in  Fig.  (5.15). 

For  example,  if  we  assume  the  maximum  field  =  45  T  at  Zw  =  0.5  m,  then  by  Eq. 
(5.13)  the  field  profile  is  given  by  Fig.  (5.14)  assuming  a  bore  diameter  of  32  mm. 

By  including  Eq.  (5.13)  in  the  simulation  we  have  the  following  results  shown  in  Fig.  5.16. 
Notice  the  results  indicate  a  final  energy  of  1.15  GeV  versus  1.2  GeV  for  a  the  static  case. 
The  top  left  graph  plots  electron  energy  gained  versus  distance  traveled.  The  bottom  left 
graph  plots  versus  distance  traveled.  The  bottom  right  graph  plots  1  —  versus  distance 
traveled.  The  only  major  difference  with  Fig.  5.2  is  that  the  electron’s  radial  position  flared 
out  in  the  last  few  millimeters.  Comparing  this  with  Fig.  5.15,  the  magnetic  field  strength 
drops  to  nearly  half  its  maximum  value.  Therefore,  the  electron  spirals  out  more. 
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Figure  5.14:  Magnetic  field  variation  along  z  with  maximum  field  af  Zw  =  0.5  m  [25] 


Figure  5.15:  Diagram  defines  di  and  62  for  calculating  the  magnetic  field  along  the  longitudinal 
axis  of  a  solenoid  represented  by  the  box. 


5.5  Diffraction  Effects 

The  numerical  solutions  presented  in  the  previous  sections  have  assumed  plane- wave  de¬ 
scription  for  the  laser  source.  In  real  lasers,  the  optical  beam  wavefronts  are  spherical  with 
the  fundamental  mode  taking  on  a  Gaussian  irradiance  pattern  transverse  to  the  optical  axis 
.  fields  for  the  electric  field  of  the  laser  where  the  wavefronts  ar.  Real  lasers  have  a  Gaussian 
approximation  to  the  laser  However,  let  us  consider  the  effects  of  Gaussian  mode.  Begin- 


67 


> 

(D 

CD 


z(m) 


Figure  5.16:  Real  magnetic  field  effects  included  in  simulation  for  Yo  =  119.  The  top  left  graph 
plots  electron  energy  gained  versus  distance  traveled.  The  top  right  graph  plots  the 
radial  position  of  the  electron  from  optical  axis.  The  bottom  left  graph  plots  versus 
distance  traveled.  The  bottom  right  graph  plots  1  —  versus  distance  traveled. 


ning  with  a  circularly-polarized  Gaussian  laser  field  in  the  fundamental  mode,  we  have  the 
following  electric  and  magnetic  field  components: 


°w(z) 

''  c  w{z) 


1 
i 

—i 

1 


(5.14) 


68 


with 


.  _ 

2R 

-”"Czr) 

(5.15) 

W  =  Wo]J 

\  {z-Zy,Y 

(5.16) 

nw^ 

ry  0 

A 

(5.17) 

R=^—{{z-z^f  +  zf^ 

Z  —  Zw 

(5.18) 

where  ^  =  krZ—  COrt,  (^g  is  the  Gaussian  phase  contribution,  w  is  the  beam  radius  radius, 
Wo  is  the  beam  waist,  is  the  Rayleigh  length,  R  is  the  phase  front  radius  of  curvature,  and 
Zw  is  the  position  of  the  laser  beam  waist. 

Substituting  Eqs.  (5.14)  into  (3.4)  then  substituting  the  real  part  into  Eq.  (3.3)  yields  the 
following: 

at  me  w 

:77  = - e  '  cos  <^G+C  l-^z) - (5-20) 

at  me  w  m 

i(yPy)-  ^^<'-'"''“%m(^.c  +  C)(l-;5=)  +  — /S,  (5.21) 

at  me  w  m 

4(rft)  =  lftcos(fc  +  C)  - ;SySin(fc  +  C)1  (5.22) 

at  me  w 

Eqs.  (5.19)  through  (5.22)  are  similar  to  Eqs.  (3.19)  through  (3.22)  in  the  previous  section 
but  they  are  now  modified  to  include  Gaussian  beam  effects  in  the  amplitude  and  phase 
terms.  To  first  order,  let  us  assume  the  following:  the  electron’s  trajectory  is  well  inside 

the  optical  mode  and  the  Rayleigh  length,  Zy,  is  two  meters  long  (twice  the  interaction 

111 

region  distance).  With  these  assumptions,  the  amplitude  factor,  Woe~''  /w,  does  not 
significantly  change  since  the  electron  radius  is  well  inside  the  optical  mode.  The  Guoy 
phase  shift,  (j)G,  is 

0G  =  -tan^^ 
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Figure  5.17:  Guoy  phase  shift  over  the  interaction  region  with  Zw  =  0.5. 


and  plotted  in  Fig.  (5.17).  As  we  expect,  to  first  order,  the  d^c/dz.  ~  about  the 

waist. 

Once  again  we  can  define  the  constant  of  the  motion,  e  =  7(1  —  based  on  Eqs.  (5.19) 
and  (5.22)  by  noting  that  the  right-hand  sides  are  identical.  Now  we  must  consider  how  (j)G 
affects  the  other  constant  of  the  motion,  rj . 

Note  the  following: 

^cos((^G+C)  =  -  sin(0G  +  C)^  +  «r(l  -/3z)sin((^G  +  C)  (5-24) 

dt  dt 

^sin(0G+C)=  cos(0G+C)^-«r(l-/3z)cos(0G  +  C)  (5-25) 

sm{(j)G  +  (r/^x)  =  ^  (r^xSin(<|)G  +  Q)  -  r/3;,cos(<^)G  +  0^  +  r^xWr(l  - /3Jcos(<|)g  +  Q 

(5.26) 

COs{(j)G+C)j^irPy)  =  ^  irPyCOs{(l)G+  Q)  +  fPy  Sm{(l)G  +  -  Y^yM^  -  l^z)  sin(<^G+  C) 

(5.27) 
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Solving  for  and  Py  in  Eqs.  (5.20)  and  (5.21)  and  substituting  in  Eq.  (5.19)  we  have 


dy 

dt 


1  /I  o 

e  B,,  dt 


■r 


d(j)G  m  dy  Eq  d 


dt  eBu  dt 


cBu  dt 


{ypyCos{(l)G  +  C))- 


Eg  d 
cBu  dt 


{ypxsm{^G  +  0)  ■ 

(5.28) 


To  first  order,  d(f>G/dt  =  {d^G/dz)  {dz/dt)  ^  —cjZr.  Substituting  this  into  Eq.  (5.28)  and 
extracting  an  overall  time  derivative  we  arrive  at 


eBu 

ygM 


(Ori^-Pzo)  +  (/3xSin(<^G  +  C)  +^yCos(<^G+  0) 

yotnc 


c  y 
IZryo 


=  constant 
(5.29) 


and  as  before  Am  =  eBu/yom  —  C0r{l  —  Pzo),  yo  =  7(0) >  £  =  7o(l  —  and  coe  =  eEolymc. 
Comparing  Eq.  (5.29)  to  Eq.  (3.44)  we  see  that  the  Gaussian  beam  description  has  added 
(j)G  in  the  trigonometric  functions  and  a  new  term,  — c7/2Z^7o.  The  magnitude  of  this  first 
order  correction  is  on  the  order  of  10^  s“^  This  is  approximately  two-orders  of  magnitude 
smaller  than  other  terms.  The  conclusion  is  that,  to  first  order,  a  Gaussian  beam  will  have 
similar  solutions  to  the  plane-wave  solutions  calculated  in  the  previous  sections. 
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CHAPTER  6: 
Conclusions 


The  relativistic  Lorentz  force  law  for  the  laser  accelerator  has  been  solved  numerically 
using  Euler-Cromer  integration  methods  for  the  first  time.  Simulation  results  indicate  that 
electron  energy  can  grow  from  about  60  MeV  to  over  1.2  GeV  in  one  meter  using  today’s 
high  powered  lasers  and  magnetic  fields.  The  electron’s  perpendicular  velocity  is  a  key 
parameter  in  achieving  good  energy  gain.  For  the  first  time,  the  complete  evolution  of 
the  electron’s  perpendicular  velocity  over  the  interaction  region  is  shown.  The  resonance 
condition  was  reexamined  and  a  minimum  energy  requirement  was  determined  for  the  first 
time.  Laser  accelerator  performance  was  examined  when  at  Am  ^  0  and  was  determined 
to  be  more  tolerant  of  off-resonance  injection  than  previously  thought. 

The  initial  electron  phase,  or  the  angle  between  J3j_  and  Er,  was  examined  and  surprisingly 
found  to  be  tolerant  of  unfavorable  initial  phase  orientation.  If  operating  near  the  minimum 
energy  and  at  resonance,  it  was  found  that  there  is  significant  energy  gain  for  all  initial 
electron  phases.  Electrons  that  were  once  thought  to  not  gain  energy  due  to  unfavorable 
initial  phase  orientations  may  still  contribute  to  a  high  energy  beam  if  the  initial  energy  is 
near  the  minimum  energy.  Therefore,  a  large  parameter  space  in  phase  and  injection  angles 
exists  for  the  laser  accelerator. 

Diffraction,  magnetic  field  solenoid  edge,  and  Guoy  phase  evolution  effects  were  examined 
and  differences  from  ideal  plane  wave  results  are  provided.  Magnetic  field  edge  effects 
were  examined  and  determined  to  not  adversely  affect  electron  energy  gain. 

A  survey  of  existing  high  powered  laser  sources,  magnetic  field  systems,  and  injectors  was 
conducted  to  give  a  sense  of  what  a  university  sized  laser  accelerator  would  have  to  be. 
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APPENDIX:  Simulation  Codes 


There  were  many  individual  programs  that  were  created  in  C  or  Matlab  for  a  specific  task 
such  as  calculating  resonance  angle  for  a  given  magnetic  field,  laser  wavelength,  and  initial 
electron  injection  energy.  In  this  appendix  a  sample  of  the  laser  accelerator  engine  code 
is  included  for  reference.  The  numerical  integration  was  programmed  in  C  with  graphics 
exported  for  use  in  Matlab. 

A.l  Laser  Accelerator  Simulation  Core 

// 

//  main.c 
//  LAsimV 
// 

//  Created  by  Ricardo  Vigil  on  3/4/14. 

//  Copyright  (c)  2014  Ricardo  Vigil.  All  rights  reserved. 

// 

#include  <stdio.h> 

#include  <stdlib.h> 

#include  <math.h> 

#include  <time.h> 

#define  pi  3.1415926535897932 

//double  pi=3 .141592653589793238; 

double  c  =  2 . 9 97 92458E8 ; 
double  uo  =  12 . 5 6637 0 614E-7 ; 
double  eo  =  8 . 854187817E-12; 
double  q  =  1 . 602176565E-19; 
charge  (C) 

double  m  =  9 . 1 0 9382 91E-31 ; 


//  speed  of  light  (m/s) 

//  Permeability  of  vacuum  (N  A^-2) 
II  permittivy  constant  (E/m) 

/ /  absolute  value  of  electron 

//  mass  of  electron  (kg) 
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double  mc2_ev  =  0 . 510998928E6; 

double  Zr  =  2.0; 

double  lambda_r  =  IE- 6; 

laser 

double  k; 

double  P  =  IE 15; 

is  lO-'lS 

double  Ipk  =  1.0E25; 

Vulcan  Laser 

double  delta_z  =  lE-5; 

size  (m) 

double  L  =  1.0; 

double  wo; 

m  for 

double  xo  =  0.0; 
double  yo  =  0.0; 
double  dt  =  7E-14; 
double  tiempo; 
double  Bu  =  45.0; 
double  width; 
double  wr; 
double  epsilon; 
double  dw; 

double  gamma_min,  gamma_o,  ggamma; 
double  alpha; 

double  A,  Eo,  beta_o,  we,  bzo,  bpo 
g_star; 

double  zeta,  bp,  bz,  bx,  by,  bz,  x 
g_check_bxby; 

double  ax,  ay,  az,  dx,  dy,  dz,  t, 
double  bxold,  byold,  bzold,  LP; 
double  phi_e,  phi_l,  bxl,  byl.  Ex, 


//  elctron  rest  mass  (eV) 

/ /  Rayleigh  length  (m) 

/ /  laser  wavelength  (m)  Vulcan 

//  wave  number  (m^-1) 

/ /  laser  peak  power  (W)  source 

//  peak  intensity  W  m"'-2  at 

//  propagation  length  step 

/ /  length  of  interest  (m) 

//  laser  beam  waist  radius  in 

//  initial  electron  position 
//  initial  electron  position 
//  time  step 

//  time  for  light  to  travel 
//  External  B  field 
/ /  pipe  width  (m) 


wco,  g_m,  bpmax,  rot_angle, 
Yr  'z-r  g_check_bp, 

,  a, b, pi , p2 , bxo,  byo; 

Ey,  brk,  Bu,  r_approx. 


76 


alpha_ 

double 

dalpha 

alpha_ 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 


o,  phi,  eta; 

gamma_top,  gamma_bot,  dgamma,  alpha_top,  alpha_bot, 

,  epsilon_res, 
res,  kappa; 

pO,  px,  py,  pz,  ggammaold,  Eng,  eta,  bp_bxby,  bp_epsilon; 
pOold,  pxold,  pyold,  pzold; 
ctheta,  stheta,  stheta2,  brk,  brk2; 
dl [1000001]  ; 

delta_g,  g,  IntG,  epsilon_f; 
frac_a,  frac_b,  frac_ep,  frac_eta; 
ao,  bo,  epsilon_o,  eta_o,  dwo; 

Ex,  Ey; 

Eangle,  bpangle; 
etap,  etam; 
theta; 

etaprime,  temp,  wtf; 
bl,  b2; 

ggammal,  zl,  test; 
aa,  bb,  cc,  zz; 
zo; 

frac_dw,  dw_anal,  dw_te,  alpha2,  gp,  gm; 
denom,  denom_ul,  denom_ul_explicit ; 
factorl,  factor2,  factorl2,  denomEactor; 
factors,  factorSeta; 
f actor4 ; 

expandedEact orS ; 
expandedDenom; 
expandedDenom_dwwe ; 
a  a_u ,  bb_u ,  c  c_u ; 

d,  f; 

dwwe_zero,  dwwe_l,  dwwe_2,  LHS,  RHS; 
delta,  rho; 
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double  deltaomega; 
double  z_anal; 


int  i,j; 
int  nalpha=150; 
int  ngamma=lll; 
int  nt; 

int  count=l 000001; 


int  flag  =  4; 


// 

case 

I . 

resonance  & 

const s 

// 

case 

2  . 

resonance  & 

no 

const s 

// 

case 

3. 

no  resonance 

& 

const s 

// 

case 

4  . 

no  resonace 

& 

no  consts 

// 

case 

5  . 

explict  use 

of 

cos (theta 

int  main (void) 

{ 

clock_t  begin,  end; 
double  time_spent; 
begin  =  clock  ( ) ; 

FILE  *outl0,  *outll; 
FILE  * out 30,  * out 31; 
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// - 

//  Initialize  arrays 

// - 

for (i=0; i<1000001; i++)  dl[i]  =  0.0; 


// - 

//  Calculate  Global  constants 

// - 

nt  =  round(L/ (c*dt) )  +  1; 

//  number  of  time  steps 
width  =  32.0  /  1000; 

//  width  of  pipe. 

wo  =  sqrt (Zr*lambda_r/pi) ; 

//  mode  radius  at  waist 
A=pi*wo*wo; 

//  area  of  laser  beam  (m"'2)  at  waist 
Eo  =  sqrt (  (P) / (c*eo*A) ) ; 

//  Electric  Eield  (V/m) 

wr  =  (2*pi*c/lambda_r) ; 

//  laser  field  angular  freq 
k  =  (2*pi) /lambda_r; 

/ /  laser  wave  number 
deltaomega  =  0.0; 

printf ( "wr\t \t\t  =  % . 1 6e\n" , wr ) ; 

printf ( "deltaomegaXt  =  %. 1 6e\n" , deltaomega) ; 

wr  =  wr  +  deltaomega/2 . 0 ; 

printf ("wr  (new) \t  =  % . 1 6e\n" , wr) ; 


// 
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//  Initial  Conditions  (Global) 

// - 

printfC'Eo  =  %  .  1 6e\n"  ,  Eo)  ; 

X  =  0.0; 

y  =  0.0; 
z  =  0.0; 
z_anal  =  0.0; 
r  =  x*x  +  y*y; 
r_approx  =  0.0; 
t  =  0.0; 

LP  =  0.0; 
rot_angle  =  0.0; 
zeta  =  k*z  -  wr*t; 

switch  (flag) { 

case  4  : 

{ 

printfC'case  4.\n\n"); 
printf ( "wo\t \t  =  % . 1 6e\n" , wo) ; 

EILE  *out40,  *out41,  *out42,  *out43,  *out44; 

Bu  =  45.0; 

out40  =  f open ( " /Users/vigilr/Dropbox/EEL/ 
LAsimV /lasim7_case4_Io . txt " , "w") ; 

out41  =  f open (" /Users/vigilr/Dropbox/EEL/ 
LAsimV /Iasim7_case4 . out ", "w" )  ; 

out42  =  f open (" /Users/vigilr/Dropbox/EEL/ 
LAsimV / Iasim7_case4_int . out" , "w" ) ; 

out43  =  f open (" /Users/vigilr/Dropbox/EEL/ 
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LAsimV /Iasim7_case4_test . out ", "w" ) ; 

out44  =  f open ( " /Users/vigilr/Dropbox/FEL/ 
LAsimV /Iasim7_case4_f racDW . out " , "w") ; 


// - 

//  Initial  Conditions  (specific  to  Case  4) 
//No  constraints  on  gamma_o,  alpha, 
and  phi  for  a  fixed  Bu . 

// - 

ggamma  =  150.0; 
ggammaold  =  ggamma; 
gamma_o  =  ggamma; 

phi  =  pi;  //3 . 1292756331662902e-01; 
alpha_res  =  acos ( (gamma_o  -  ( (q*Bu) / (wr*m) )  ) / 

(sqrt (gamma_o*gamma_o  -  1) )  ) ; 

alpha  =  alpha_res; 


bz  =  1 . 0/ggamma*sqrt (ggamma*ggamma  -  1 . 0 ) *cos (alpha) ; 
printf ( "bzo\t \t=  % . 1 6e\n" , bz ) ; 

bp  =  1 . 0/ggamma*sqrt (ggamma*ggamma  -  1 . 0 ) *sin  (alpha) ; 
printf ( "bpo\t\t=  % . 1 6e\n" , bp) ; 

bx  =  bp*cos  (phi) ;  //  eq  4 

by  =  bp*sin(phi);  //  eq  5 

beta_o  =  1 . 0/ggamma*sqrt (ggamma*ggamma  -  1.0); 
if  (fabs (bp)  <  lE-15)  bp=0.0; 
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if  (fabs  (bx)  <  lE-15)  bx=0.0; 
if  (fabs (by)  <  lE-15)  by=0.0; 
zo  =  atari  (bx/by) /k; 

epsilon  =  ggamma* ( 1-bz ) ;  //  eg  6 


bp_bxby  =  sqrt (bx*bx  +  by*by) ; 
bp_epsilon  =  (1/ggamma)* 
sqrt (2*ggamma*epsilon  -  epsilon*epsilon  -  1); 


bpo  = 

bp_epsilon; 

bxold 

=  bx; 

byold 

=  by; 

bzold 

=  bz ; 

pO  = 

ggamma*m*c; 

// 

eq 

10 

px  = 

ggamma *m*c*bx; 

// 

eq 

7 

PY  = 

ggamma *m*c*by; 

// 

eq 

8 

pz  = 

ggamma*m*c*bz ; 

// 

eq 

9 

pOold 

II 

P 

o 

pxold 

=  px; 

pyold 

II 

P 

P 

pzold 

=  pz; 

g_check_bxby  =  l/sqrt(l 

-  bx*bx  - 

>1 

b 

1 

*by 

ax  = 

c*bx/ dt ; 

ay  = 

c*by/ dt ; 

az  = 

c*bz / dt ; 

dx  = 

c*by*az  -  ay*c*bz; 

dy  = 

ax*c*bz  -  c*bx*az; 

dz  = 

c*bx*ay  -  ax*c*by; 

zeta 

=  0.0;//k*z  -  wr*t; 

Ex  = 

cos (zeta) ; 

Ey  = 

-1 . 0*sin (zeta)  ; 

fprintf (out 43 , "% . 1 6e\n" , 

(pxold*cos  (zeta 

-  bz *bz ) ; 


-  pyold  * 


// 


sin ( zeta) 
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//- 

— 

— 

— 

— 

— 

— 

— 

— 

// 

initial  angles 

//- 

— 

— 

— 

— 

— 

— 

— 

— 

if 

(Ex 

> 

0 

&& 

Ey 

> 

0  ) 

Eangle  = 

atan (Ey/Ex) ; 

if 

(Ex 

< 

0 

&& 

Ey 

> 

0  ) 

Eangle  = 

atan (Ey/Ex) 

+  pi; 

if 

(Ex 

< 

0 

&& 

Ey 

< 

0  ) 

Eangle  = 

atan (Ey/Ex) 

+  pi; 

if 

(Ex 

> 

0 

&& 

Ey 

< 

0  ) 

Eangle  = 

atan (Ey/Ex) 

+  2.0  *pi ; 

if 

(bx 

> 

0 

&& 

by 

> 

0  ) 

bpangle 

=  atan (by/bx) 

r 

if 

(bx 

< 

0 

&& 

by 

> 

0  ) 

bpangle 

=  atan (by/bx) 

+  pi; 

if 

(bx 

< 

0 

&& 

by 

< 

0  ) 

bpangle 

=  atan (by/bx) 

+  pi; 

if 

(bx 

> 

0 

&& 

by 

< 

0  ) 

bpangle 

=  atan (by/bx) 

+  2.0  *pi ; 

theta  =  Eangle  -  bpangle; 

if  (theta  <  0  )  theta  =  theta  +  2.0*pi; 

if  (theta  >  0  &&  theta  >  pi)  theta  =  2.0*pi  -  theta; 

if  (phi>0  &&  phi<pi)  theta  =  phi; 

// - 

/ /  Calculate  derived  parameters 

// - 

test  =  bx*cos  (k*z-wr*t+k*zo)  -  by*sin (k*z-wr*t+k*zo) ; 

we  =  (q*Eo) / (gamma_o*m*c) ; 

wco  =  (q*Bu) / (m*ggamma) ; 

g_star  =  (epsilon*epsilon+l ) /epsilon; 

gamma_min  =  (epsilon*epsilon+l ) /  (2 *epsilon) ; 

bpmax  =  (epsilon/sqrt (epsilon*epsilon+l ) ) 

dw  =  wco  + 

(wr/ggamma) *sqrt (ggamma*ggamma  -  1 ) *cos (alpha)  -  wr;  //  eg  23 
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dw_anal  =  (q*Bu) / (ggamma*m)  + 

wr*bz  -  wr; 

a  =  ggamma*bx  -  ( (g*Eo) / (m*c*wr) ) *sin ( zeta) 

+  ( (g*Bu) / (m*c) ) *y; 

b  =  ggamma*by  -  ( (g*Eo) / (m*c*wr) ) *cos ( zeta) 

-  ( (g*Bu) / (m*c) ) *x; 

temp  =  (bx*sin  (zeta)  +  by*cos  ( zeta) ) ; 
eta  =  ggamma*dw  +  ggamma*we* (bx*sin ( zeta) 

+  by*cos (zeta) ) ;  //  eg  24 

etap  =  ggamma*dw  + 
ggamma*we*bp_epsilon*sin (theta) ; 

etam  =  ggamma*dw  - 
ggamma*we*bp_epsilon*sin (theta)  ; 

wtf  =  asin ( (wr  -  wco  -  wr*beta_o*cos (alpha) ) 

/ (we*beta_o*sin (alpha) ) )  ; 

stheta  =  (eta  -  ggamma*dw) / 

(we*sgrt (2*ggamma*epsilon  -  epsilon*epsilon  -  1)); 

stheta2  =  (bx*sin ( zeta)  +  by*cos (zeta) ) /bp_bxby; 

ctheta  =  (bx*cos ( zeta)  -  by*sin ( zeta) ) /bp_epsilon; 

brk  =  bx*cos(zeta)  -  by*sin ( zeta) ;  //bp  *  cos (theta) ; 

brk2  =  bx*sin(zeta)  +  by*cos ( zeta) ;  //bp  *  sin (theta) ; 

f rac_a  =  0.0; 

frac_b  =  0.0; 

frac_ep  =  0.0; 

frac_eta  =  0.0; 

3.0  =  di} 

bo  =  b; 

epsilon_o  =  epsilon; 
eta_o  =  eta; 

// - 

//  Inputs 
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// - 

fprintf (out 4 0 , "case  4.  no  resonance  &  consts\n"); 
fprintf  (out40,  "\t - Inputs - \n")  ; 


// - 

//  Equations  of  Motion 
// - 


fprintf (out 4 1 ,"%.16e\t%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t%.16e\t 
% . 16e\t%. 16e\t%. 16e\t%16\n", 
t , X, y , z , r, r_approx, ggamma, bx, by, bz , 
bp_bxby , bp_epsilon, g_check_bxby , LP , zeta, 

St beta, ctheta, brk, epsilon, a, b, eta, 
frac_a, frac_b, frac_ep, frac_eta, 
brk2 , stheta2 , theta, etap, etam, 
et aprime , dw_anal , z_anal ) ; 

for (i=0; i<nt; i++) { 

II  — 

pO  =  pO  + 

( (-q*Eo/ (ggammaold*  m*c) ) * (pxold*cos ( zeta)  -  pyold  *  sin (zeta)  ) 

ggamma  =  pO/ (m*c) ; 

bz  =  (ggamma  -  epsilon) /ggamma; 

z  =  z  +  bz*c*dt; 


)  *dt; 
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// - 

px  =  pxold  +  ( (-q*Eo) * (epsilon/ggammaold) *cos ( zeta) 
-  (q*Bu*pyold) / (ggammaold*m) ) *dt ;  //  eq  15 
bx  =  px/ (ggamma*m*c) ; 

X  =  X  +  bx*c*dt; 

// - 

PY  =  pyold  +  ( (q*Eo) * (epsilon/ggammaold) *sin ( zeta) 

+  (q*Bu*pxold) / (ggammaold*m) ) *dt ;  //  eq  18 
by  =  py/ (ggamma*m*c) ; 
y  =  y  +  by*c*dt; 

// - 

pxold  =  px; 
pyold  =  py; 
ggammaold  =  ggamma; 

// - 

t  =  t  +  dt ; 

zeta  =  k*z  -  wr*t; 

// - 

bp_bxby  =  sgrt (bx*bx  +  by*by) ; 

bp_epsilon  =  ( 1 /ggamma) *sqrt (2 *ggamma*epsilon 
-  epsilon*epsilon-l ) ; 

if  ( (2*ggamma*epsilon  -  epsilon*epsilon-l )  <  0  )  { 

bp_epsilon  =  0.0; 

} 

// - 

stheta  =  (eta_o  -  ggamma*dw) / 

(we*sqrt (2*ggamma*epsilon  -  epsilon*epsilon  -  1));; 

stheta2  =  (bx*sin ( zeta)  +  by*cos ( zeta) ) /bp_bxby; 
ctheta  =  (bx*cos ( zeta)  -  by*sin ( zeta) ) /bp_epsilon; 
brk  =  bx*cos(zeta)  -  by*sin  ( zeta) ;  // 
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brk2  =  bx*sin(zeta)  +  by*cos (zeta)  ; 
epsilon_f  =  ggamma* ( 1-bz ) ; 


//- 

— 

— 

- 

— 

— 

— 

// 

Angles 

//- 

— 

— 

- 

— 

— 

— 

Ex 

=  < 

SOS 

( 

zeta) ; 

Ey 

=  -i 

sin  ( 

zeta) ; 

if 

(Ex 

> 

0 

&& 

Ey 

>  0 

)  Eangle  = 

atan (Ey/Ex)  ; 

if 

(Ex 

< 

0 

&& 

Ey 

>  0 

)  Eangle  = 

atan (Ey/Ex) 

if 

(Ex 

< 

0 

&& 

Ey 

<  0 

)  Eangle  = 

atan (Ey/Ex) 

if 

(Ex 

> 

0 

&& 

Ey 

<  0 

)  Eangle  = 

atan (Ey/Ex) 

if 

(bx 

> 

0 

&& 

by 

>  0 

)  bpangle 

=  atan  (by/bx) 

if 

(bx 

< 

0 

&& 

by 

>  0 

)  bpangle 

=  atan (by/bx) 

if 

(bx 

< 

0 

&& 

by 

<  0 

)  bpangle 

=  atan (by/bx) 

if 

(bx 

> 

0 

&& 

by 

<  0 

)  bpangle 

=  atan  (by/bx) 

theta  =  Eangle  -  bpangle; 

if  (theta  <  0  )  theta  =  theta  +  2.0*pi; 

if  (theta  >  0  &&  theta  >  pi)  theta  =  2.0*pi 

bl  =  bp_epsilon*cos (theta)  ; 
b2  =  bp_epsilon*sin (theta) ; 

II  — 

a  =  ggamma*bx  -  ( (q*Eo) / (m*c*wr) ) *sin  (zeta) 

+  ( (q*Bu) / (m*c) ) *y; 


pi; 

pi; 

2 . 0  *pi ; 

+  pi; 

+  pi; 

+  2.0  *pi 


theta; 
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b  =  ggamma*by  -  ( (q*Eo) / (m*c*wr) ) *cos (zeta) 

-  ( (q*Bu) / (m*c) ) *x; 

eta  =  ggamma* (dw  +  ( (q*Eo) / (gamma_o*m*c) ) * 

(bx*sin ( zeta)  +  by*cos ( zeta) ) )  ; 

etap  =  ggamma* (dw  +  ( (q*Eo) / (gamma_o*m*c) ) * 
bp_epsilon*sin (theta) )  ; 

etam  =  ggamma* (dw  -  ( (q*Eo) / (gamma_o*m*c) ) * 
bp_epsilon*sin (theta) )  ; 

etaprime  =  ggamma  +  ggamma* (Eo/c*Bu) * 

(bx*sin  ( zeta)  +  by*cos ( zeta) )  -  (m*epsilon*wr*ggamma) / (q*Bu) 

frac_a  =  (a  -  ao)/ao; 

frac_b  =  (b  -  bo) /bo; 

frac_eta  =  (eta  -  eta_o) /eta_o; 

frac_ep  =  (epsilon_f  -  epsilon) /epsilon; 

II  — 

r  =  sqrt (x*x  +  y*y)  ; 

r_approx  =  (m*ggamma*sqrt (bx*bx+by*by ) *c) / (q*Bu) ; 

// - 

ax  =  c*  (bx  -  bxold)/dt; 
ay  =  c* (by  -  byold)/dt; 
az  =  c*  (bz  -  bzold)/dt; 
bxold  =  bx; 
byold  =  by; 
bzold  =  bz; 

dx  =  c*by*az  -  ay*c*bz; 

dy  =  ax*c*bz  -  c*bx*az; 

dz  =  c*bx*ay  -  ax*c*by; 

LP  =  ( (uo*q*q*ggamma*ggamma*ggamma*ggamma* 

ggamma* ggamma) / (6*pi*c)  ) * 

(ax*ax  +  ay*ay  +  az*az  - 

(1/  (c*c) ) * (dx*dx  +  dy*dy  +  dz*dz) ) ; 
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// - 

rot_angle  =  rot_angle  +  (q*Bu) / (m*ggamma) *dt ; 

// - 

g_check_bxby  =  l/sqrt(l  -  bx*bx  -  by*by  -  bz*bz) ; 

// - 

dw_anal  =  (q*Bu) / (ggamma*m)  -  wr*(l-bz); 

// - 

z_anal  =  (c/gamma_o)* 

(1/  (3*we*we*we*we*epsilon*epsilon) ) * 
sqrt (2*we*we*epsilon*ggamma  - 
we*we*epsilon*epsilon  -  we*we  -  eta*eta) * 
(we*we*epsilon*ggamma  - 
we*we*epsilon*epsilon  +  we*we  + 
eta*eta  -  3*we*we*epsilon*epsilon)  ; 


// - 

//  Matlab  output 
// - 


fprintf (out41, "%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%  .  16e\t%  .  16e\n’’, 
t , X , y , z , r , r_appr ox , ggamma , 
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bx , by , b  z , bp_bxby , 

bp_epsilon, g_check_bxby , LP , 

zeta, St beta, ctheta, brk, epsilon_f , 

a, b, eta, f rac_a, f rac_b, f rac_ep, f rac_eta, 

brk2 , sthet a2 , theta, etap, et am, etaprime , dw_anal , z_anal ) ; 
} 

// - 

/ /  Outputs 

// - 

epsilon_f  =  ggamma* ( 1-bz ) ; 
a  =  ggamma*bx  - 

( (q*Eo) / (m*c*wr) ) *sin (zeta)  +  ( (q*Bu) / (m*c) ) *y; 

b  =  ggamma*by  - 

( (q*Eo) / (m*c*wr) ) *cos  (zeta)  -  (  (q*Bu) / (m*c) ) *x; 

eta  =  ggamma* (dw  + 

( (q*Eo) / (ggamma*m*c) ) * (bx*sin (zeta)  +  by*cos (zeta) ) ) ; 


fprintf (out40, "\n\t - Calculated  Values - \n") ; 

fprintf (out40, "?_max\t= 

% . 1 6e\n?pmax\t=  % . 1 6e\n" , g_m, bpmax) ; 
fprintf  (out40, 

"\n\t - Outputs - \n")  ; 

fprintf (out 4 0 ," ?_E\t=  %.16e\n 
Eng_E\t=  %.16e\n?z\t=  %.16e\n 
?p\t=  % . 1 6e\nRev\t=  %.16e\n 
r_E\t=  % . 1 6e\nr_apprx\t=  %.16e\n 
L\t=  % . 1 6e\nTime\t=  %.16e\n", 
ggamma, mc2_ev*ggamma/ 1E9 , bz , 
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bp, rot_angle/ (2*pi) , r, r_approx, z, t) ; 

print f ( " ?_star\t=  % . 1 6e\n?_min\t=% . 1 6e\n 
?_o\t\t=  % . 1 6e\n?_F\t \t=  %.16e\n 
Eng_F\t=  % . 1 6e\n? \t \t=  %.16e\n 
?\t\t=  % . 1 6e\nBu\t \t=  %.16e\n 
?\t\t=  % . 1 6e\na\t \t=  %.16e\n 
b\t\t=  % . 16e\n?\t\t=  %.16e\n 
we\t\t=  % . 1 6e\nwc\t \t=  %.16e\n 
wr\t\t=  % . 1 6e\n?w\t \t=  %.16e\n 
Revs\t=  %g\nnt\t\t=  %d\nz\t\t=  %.16e\n", 
g_star, gamma_min, gamma_o, 
ggamma, mc2_ev*ggamma/ 1E9 , alpha, 
phi, Bu, epsilon, a, b, eta_o, 
we, wco, wr, dw, rot_angle/ (2*pi) ,nt,z) ; 
fclose (out40) ; 
fclose (out41) ; 
fclose  (out42 ) ; 
fclose  (out43) ; 
fclose (out44 ) ; 
break; 

} 

default : 

printfC'try  again\n"); 

}  / /  end  switch 

// - 

end  =  clock  ( ) ; 

time_spent  =  (double)  (end  -  begin)  /  CLOCKS_PER_SEC; 
printf ( " \n\nElapsed :  %.16g  secondsXn",  time_spent) ; 

return  0; 
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} 


A.2  Laser  Accelerator  Matlab  Graphics  Script 

clear 
close  all 
clc 

load  Iasim7_case4 . out ; 
load  Iasim7_case4_int . out ; 
load  Iasim7_case4_l Om . out ; 

c  =  2 . 99792458E8;  % 

uo  =  12 . 566370614E-7;  % 

eo  =  8 . 854187817E-12;  % 

q  =  1 . 602176565E-19;  % 

m  =  9 . 10938291E-31;  % 

mc2_ev  =  0 . 510998928E6;  % 

Eo  =  4 . 3401054911612872e+ll; 
we  =  (q*Eo) / (119 . 0*m*c) ; 


speed  of  light  (m/s) 

Permeability  of  vacuum  (N  A^-2) 
permittivy  constant  (E/m) 
absolute  value  of  electron  charge  (C) 
mass  of  electron  (kg) 
elctron  rest  mass  (eV) 

%  1/s 


%zlO  =  Iasim7_case4_l Om ( : , 4 ) ; 
%ggammalO  =  Iasim7_case4_l Om ( : , 7 ) ; 


t  =  Iasim7_case4 ( : ,  1 )  ; 

X  =  Iasim7_case4 ( : ,  2 )  ; 
y  =  Iasim7_case4 ( : ,  3 )  ; 
z  =  Iasim7_case4 ( : ,  4 )  ; 
r  =  Iasim7_case4 ( : ,  5 )  ; 
r_a=  Iasim7_case4 ( :  ,  6 )  ; 
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ggamma  =  Iasim7_case4 ( : , 7 ) ; 
bx  =  Iasim7_case4 ( : , 8 ) ; 
by  =  Iasim7_case4 ( : , 9 ) ; 
bz  =  Iasim7_case4 ( : , 1 0 ) ; 
bp_bxby  =  Iasim7_case4 ( : , 1 1 )  ; 
bp_epsilon  =  Iasim7_case4 ( : , 12 ) ; 
g_check_bxby  =  Iasim7_case4 ( : , 13) ; 
LP  =  Iasim7_case4 ( :  ,  14 )  ; 
zeta  =  Iasim7_case4 ( :  ,  15)  ; 
stheta  =  Iasim7_case4 ( :  ,  1 6 )  ; 
ctheta  =  Iasim7_case4 ( :  ,  17 )  ; 
brk  =  Iasim7_case4 ( : , 1 8 ) ; 
epsilon  =  Iasim7_case4 ( : , 1 9)  ; 
a  =  Iasim7_case4 ( : ,  2 0 )  ; 
b  =  Iasim7_case4 ( : ,  2 1 )  ; 
eta  =  Iasim7_case4 ( :  ,  22 )  ; 
frac_a  =  Iasim7_case4 ( : , 23)  ; 
frac_b  =  Iasim7_case4 ( : , 24 ) ; 
frac_ep  =  Iasim7_case4 ( : , 25) ; 
frac_eta  =  Iasim7_case4 ( : , 2 6 ) ; 
brk2  =  Iasim7_case4 ( : , 27 ) ; 
stheta2  =  Iasim7_case4 ( :  ,  2 8 )  ; 
theta  =  Iasim7_case4 ( : , 2 9 )  ; 
etap  =  Iasim7_case4 ( :  ,  30 )  ; 
etam  =  Iasim7_case4 ( : , 31 ) ; 
etaprime  =  Iasim7_case4 ( : , 32 ) ; 
dw_anal  =  Iasim7_case4 ( : , 33) ; 
z_anal  =  Iasim7_case4 ( :  ,  34 )  ; 

D  =  z_anal  -  ( z_anal (2 ) -z_anal ( 1 ) ) ; 
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zl  = 

lasim7_ 

_case4 

_int  (  : 

gl  = 

lasim7_ 

_case4 

_int  (  : 

ao  = 

a 

(1)  ; 

bo  = 

b 

(1) ; 

f  rac_ 

_a 

=  (a 

-  ao) 

/ao; 

f  rac_ 

_b 

=  (b 

-  bo) 

/bo; 

figure  ( 1 ) 

subplot  (3,2,1) 

set (gca, 'Font Size', 24) 

plot (z, ggamma*mc2_ev/lE9, ' LineWidth' 

%hold  on 

%plot ( z , g_check_bxby*mc2_ev/ 1E9, ' r' , 

xlabel (' z' , 'EontSize' ,24) 

ylabel  ( '  \gammamc'^2  (GeV)  '  ) 

xlim  (  [ 0  max ( z ) ] ) 

grid 

legend (' gamma' ,' gamma  check') 

subplot (3,2,2) 

%plot ( z , r, ' LineWidth'  ,  3 ) 
set (gca, 'EontSize', 24) 

%hold  on 

plot ( z , r_a, ' r ' , ' LineWidth' , 2 ) 
xlabel ( ' z  (m) ' ) 
ylabel ( ' r  (m)  '  ) 
grid 

xlim ( [0  max ( z ) ] ) 


,3) 

'  LineWidth' , 2 ) 
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legend ( '  r 


r  approx' ) 


subplot  (3,2,3) 

set (gca, 'Font Size', 24) 

plot ( z , bp_epsilon, ' LineWidth' , 3 ) 

%hold  on 

%plot ( z , bp_bxby , ' r' , ' LineWidth'  ,  2 ) 

%hold  on 

%plot ( z , brk, ' g' ) 

xlabel  ( ' z ' ) 

ylabel (' \beta\perp' ) 

xlim ( [0  max ( z ) ] ) 

grid 

subplot  (3,2,4) 

set (gca, 'Font Size', 24) 

plot ( z , 1-bz , ' LineWidth' , 3 ) 

xlabel  ( ' z ' ) 

ylabel ('1  -  \beta_z') 

grid 

xlim ( [0  max ( z ) ] ) 

subplot  (3, 2, 5 : 6) 

plot ( z , s theta, ' LineWidth' , 3) 

hold  on 

plot ( z , ctheta, ' r ' , ' LineWidth' , 3 ) 

set (gca, 'Font Size', 24) 

xlabel  ( ' z ' ) 

ylabel ( ' sin/cos ' ) 

xlim([0  1]) 

legend  ('  sin  (theta)  '  ,  'cos  (theta)  '  ,  'brk'  ) 
title ('does  \theta  =  \theta_o' ) 
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figure  ( 2 ) 

set (gca, 'Font Size' ,24) 

plots  (x, y, z ) 

xlabel  ( ' x' ) 

ylabel ( ' y' ) 

zlabel  ( ' z ' ) 

grid 


figure  ( 3 ) 

subplot  (2,2,1) 

plot (z, frac_ep) 

set (gca, 'Font Size' ,24) 

xlabel  ( ' z  (m) ' ) 

ylabel ('frac  ep' ) 

xlim ( [0  max ( z ) ] ) 

subplot  (2,2,2) 

plot ( z , f rac_a) 

set (gca, 'Font Size', 24) 

xlabel  ( ' z  (m) ' ) 

ylabel  ( ' frac_a' ) 

xlim ( [0  max ( z ) ] ) 

subplot (2,2,3) 

plot ( z , f rac_b) 

set (gca, 'Font Size', 24) 

xlabel  ( ' z  (m) ' ) 

ylabel ( ' frac_b' ) 

xlim ( [0  max ( z ) ] ) 

subplot (2,2,4) 

plot  (z, frac_eta) 

set  (gca, 'Font Size'  ,24) 
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xlabel ( ' z  (m) ' ) 
ylabel('frac  eta') 
xlim ( [0  max ( z ) ] ) 

figure ( 4 ) 

plot ( z , brk, z , brk2 ) 
legend ( ' brk' , ' brk2 ' ) 
xlim ( [0  max ( z ) ] ) 

figure  ( 5 ) 

plot (z, stheta, z, stheta2) 
legend ( ' sthetal ' , ' stheta2 ' ) 
xlim ( [0  max ( z ) ] ) 

figure ( 6 ) 
plot  ( z , theta) 
xlim  (  [0  0.1]) 
ylim  (  [0  2*pi ] ) 

figure  ( 7 ) 

plot ( z , eta, z , etap, z , etam) 
legend ( ' eta' , ' etap' , ' etam' ) 
xlim ( [0  max ( z ) ] ) 

figure ( 8 ) 

plot  ( z , ggamma, ' LineWidth'  ,  3) 

set (gca, 'Font Size'  ,26) 

set  (gcf, 'color' ,  'w'  )  ; 

xlim ( [0  max ( z ) ] ) 

xlabel ( ' z  (m)  '  ) 

ylabel (' \gamma'  ) 

grid 
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figure  ( 9 ) 

set (gcf, 'color' , 'w' ) ; 
set (gca, 'FontSize' , 32)  ; 

plot (z_anal, ggamma*mc2_ev/lE9, ' k' , ' LineWidth' , 5) 
hold  on 

plot (z, ggamma*mc2_ev/lE9, ' — r' , ' LineWidth' , 2) 
hold  on 

%plot (zlO, ggammalO*mc2_ev/lE9, ' g' , ' LineWidth' , 3) 

xlabel ( ' z  (m) ' ) 

ylabel  ( '  \gamma  mc"'2  (GeV)  '  ) 

legend (' Analytic' , ' Simulation' ) 

xlim ( [0  max ( z ) ] ) 

ylim ( [min (ggamma) *mc2_ev/ 1E9*  0 . 7  1 . l*max (ggamma) *mc2_ev/ 1E9 ] ) 
grid 
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